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FIRST  BATCH  OF  CHINESE  PATENTS  TO  BE  ANNOUNCED  IN  SEPTEMBER 

Tianjin  JISHU  SHICHANG  BAO  [TECHNOLOGY  MARKET  WEEKLY]  in  Chinese  3  Sep  85 

p  1 

[Article  by  special-assignment  reporter  Shen  Jiajian  [3947  0857  1670]:  "The 
First  Batch  of  Chinese  Patents  Will  be  Announced  Within  the  Month;  After 
the  Inventions  are  Announced  They  Will  Enter  the  Objection  Stage  and  Receive 
Temporary  Legal  Protection,  with  Uncompensated  Use  Not  Permitted." 

[Text]  Director  of  the  Chinese  Patent  Bureau  Huang  Kunyi  [7806  0981  4135] 
announced  at  the  recently  assembled  first  national  patent  documentation 
conference  that  the  first  batch  of  Chinese  patents  would  be  announced  in 
September,  to  let  the  public  carry  out  discussion  (the  period  for  objections 
being  3  months),  and  that  more  patents  would  be  announced  each  month. 

The  implementation  of  our  country's  patent  law  will  enter  a  new  stage. 

Huang  Kunyi  said  the  major  work  of  this  stage  is  to  announce  how  to  treat 
the  objection  process  and  the  openness  under  legal  protection.  The  objection 
process  is  an  examination  by  the  masses,  weighing  the  inventions  on  the 
scales  of  international  standards,  so  that  only  the  inventions  that  really 
pass  will  be  approved.  Moreover,  the  public  should  know  that  the  openness 
under  legal  protection  is  different  from  the  former  openness  without  legal 
protection.  The  temporary  legal  protection  begins  from  the  day  of  announce¬ 
ment,  so  no  announced  invention  may  be  used  without  compensation.  This 
sort  of  common  knowledge  must  be  well  proclaimed  from  the  standpoint  of 
both  law  and  morality. 

The  first  batch  of  patents  announced  will  be  approved  in  December  after 
going  through  mass  discussion.  ?  ,  j  r'-'/'L--;,  ;  ' 

13044/12899 
CS0:  4008/2001 


1 


JPRS-CST-85-040 
20  November  1985 


NATIONAL  DEVELOPMENTS 


CAS  PATENT  INSTITUTE  OPENS  TO  OUTSIDERS 

Beijing  BEIJING  KEJIBAO  [BEIJING  S&T  WEEKLY]  in  Chinese  26  Aug  85  p  1 

[Article  by  staff  reporter  Hao  Meifang  [6727  2734  3364]:  "China  Academy 
of  Sciences  Patent  Institute  Begins  External  Operations"] 

[Text]  The  China  Academy  of  Sciences  (CAS)  Patent  Institute  officially 
began  external  operations  on  16  August  at  CAS  and  at  (Kehai  Gongsi)  in 
Zhongguan  Cun  in  Haidian  District  of  Beijing. 

This  patent  institute  is  composed  of  patent  representatives  and  patent 
workers  from  each  of  the  CAS  institutes  in  the  Beijing  area.  The  great 
majority  of  them  come  from  the  first  line  of  scientific  research,  a  con¬ 
tingent  with  a  broad  coverage-  of  topics  and  specialties,  able  to  enter  the 
ranks  of  those  who  engage  in  patent  representation  in  patent  and  economic 
legal  matters. 

Since  March  of  this  year,  the  patent  institute  has  on  the  one  hand  been 
using  all  forms  of  activities  to  enter  government  offices,  research  insti¬ 
tutes  and  laboratories  to  publicize  the  meaning  of  the  patent  system  estab¬ 
lished  by  our  country,  procedures  for  applying  for  patents,  the  method  of 
selecting  documents  and  conducting  patent  searches;  they  have  also  compiled 
all  kinds  of  patent  materials  to  introduce  to  scientific  research  personnel 
trade  knowledge  related  to  patents  and  patent  permissions.  On  the  other 
hand,  they  have  actively  helped  CAS's  own  inventors  handle  applications, 
so  that  when  the  Chinese  Patent  Law  took  effect  on  1  April,  78  applications 
were  made  to  the  Patent  Bureau  through  them,  making  up  about  2  percent 
of  the  whole  country's  applications  in  the  same  period. 

At  present,  the  patent  institute  cannot  only  take  on  patent  representation 
and  patent  legal  work  for  work  units  and  individuals  in  the  academy,  but 
also  fully  possesses  all  the  requirements  for  serving  inventors  from  all 
segments  of  the  public  with  application  and  consultation  services. 

President  of  the  CAS  Lu  Jiaxi  [4151  0857  6932],  CAS  Secretary  General  Gu 
Yijian  [7357  0110  0256],  Patent  Bureau  Director  Shen  Yaozeng  [3088  1031 
2582],  officials  of  various  relevant  ministries  and  commissions  of  the 
State  Council,  and  comrades  from  the  China  Trade  Promotion  Association 
and  the  Beijing  Patent  Administration  Bureau  attended  the  opening 
ceremony. 
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SYSTEM  OF  ENGAGING  S&T  PERSONNEL  DISCUSSED 

Beijing  BEIJING  KEJIBAO  [BEIJING  S&T  WEEKLY]  in  Chinese  2  Sep  85  p  1 

[Article  by  staff  reporter  Zhi  Qiang  [1807  1730]:  "The  State  Science  and 
Technology  Commission's  Science  and  Technology  Cadre  Bureau's  Responsible 
Person  Indicated  That  Arranging  the  Implementation  of  the  Appointment  System 
for  Personnel  Must  Be  Cautiously  Undertaken"] 

[Text]  Not  long  ago,  at  the  seminar  on  the  "appointment"  work  system  for 
S&T  assembled  by  the  State  Science  and  Technology  Commission,  a  responsible 
person  from  the  Science  and  Technology  Cadre  Bureau  stressed  that  implement¬ 
ing  the  appointment  system  must  be  cautiously  undertaken  with  reference  to 
the  problems  of  others  and  must  be  based  on  the  positive  effect  of  trans¬ 
ferring  each  category  of  personnel. 

The  responsible  person  from  the  Science  and  Technology  Cadre  Bureau  of  the 
State  Science  and  Technology  Commission  said  to  use  a  lifelong  appointment 
system  for  those  few  famous  experts  with  a  reputation  for  good  character 
and  deep  academic  contributions,  who  are  well  known  at  home  and  abroad. 

For  other  experts  suited  to  a  retirement  leave,  carry  out  a  retirement  leave 
system.  Those  who  are  still  healthy  after  retirement  who  can  persist  in 
work  can,  under  conditions  they  agree  to,  continue  fulfilling  their  function, 
ensuring  that  they  will  not  only  have  something  to  rely  on  in  their  old 
age,  but  also  something  to  do  in  old  age,  and  something  to  delight  in  in 
old  age.  This  embodies  the  spirit  of  honoring  the  old  and  respecting  the 
wise,  and  is  able  to  make  even  more  middle-aged  and  young  S&T  personnel 
become  academic  leaders,  or  to  be  appointed  to  higher  positions  better 
suited  to  them.  ;  ,.,.v  ~ 

This  responsible  person  continued  by  saying  that  personnel  who  have  not 
yet  been  appointed  must  be  correctly  treated  and  their  appointments  appro¬ 
priately  arranged.  First  they  must  have  a  correct  understanding  that 
they  are  iostly  useful  talents,  not  a  burden,  and  at  the  same  time  conditions 
must  be  created  to  arrange  for  them  to  enter  posts  in  which  they  are  more 
able  to  fulfill  their  function.  Secondly,  the  purpose  of  transferring 
personnel  is  to  better  fulfill  the  function  of  personnel  and  to  adjust 
for  surpluses  and  shortages,  which  is  a  direction  of  reform  for  the 
personnel  system,  arid  the  outmoded  idea  must  be  destroyed  that  being 
transferred  to  another  post  in  the  same  or  a  different  work  unit  is  a  sign 
of  "making  a  mistake"  or  "losing  face."  According  to  our  understanding. 
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the  personnel  sent  out  from  high-level  schools,  research  institutes,  or 
hospitals  are  science  and  technology  personnel  of  above  median  ability, 
and  not  a  few  of  them  become  backbones  of  the  organization  in  their  new 
posts,  and  some  enter  the  leading  echelons,  fulfilling  an  even  greater 
function  than  in  the  original  work  unit. 

The  responsible  person  of  the  Science  and  Technology  Cadre  Bureau  of  the 
State  Science  and  Technology  Commission  also  indicated  that  for  individuals 
who  had  not  performed  well  and  who  had  not  been  appointed  and  were  difficult 
to  transfer,  there  should,  based  on  their  educational  level,,  be  assistance 
for  their  transfers  supplemented  by  necessary  administrative  and  economic 
measures.  Some  work  units  are  strict  in  regard  to  these  administrative 
and  economic  measures,  some  are  lenient.  We  think  that  since  it  is  all 
based  on  reeducation,  using  these  methods  is  necessary,  but  they  must  be 
suitable.  The  Shenyang  Metals  Institute  of  the  Chinese  Academy  of  Sciences's 
method  is  to  withhold  1  percent  of  the  salary  every  month,  which  puts  a 
certain  pressure  on  a  person,  without  intensifying  contradictions.  This 
method  is  worthy  of  study. 
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FUJIAN  SCIENCE  AND  TECHNOLOGY  DEVELOPMENT  PLAN  DISCUSSED 
Fuzhou  FUJIAN  RIBAO  in  Chinese  26  Aug  85  p  1 

[Article  by  staff  reporter  Chen  Min zhen  [7115  2404  4176]:  "To  Work  Out 
a  Provincial  Science  and  Technology  Development  Plan  with  Special  Character is 
tics,  Almost  100  Famous  Experts  and  Professors  Assembled  in  Fuzhou"] 

[Text]  From  21-23  August,  almost  100  famous  experts  and  professors  from 
our  province  assembled  at  the  Xihu  Hotel  in  Fuxhou  to  discuss  the  question 
of  provincial  scientific  and  technological  development  strategy. 

They  mentioned  that  our  province  is  located  in  the  subtropics,  has  the 
advantages  of  "mountains,  sea,  overseas  Chinese,  and  Special  Economic  Zones," 
and  should  work  out  a  S&T  development  plan  with  the  province's  special 
characteristics.  S&T  and  the  economy  are  inseparably  joined,  and  when  the 
S&T  development  plan  is  worked  out  it  should  be  adapted  to  economic  develop¬ 
ment.  They  hoped  to  form  leading  "locomotive"  industries  among  the  numerous 
industries,  and  suggested  placing  the  biotechnology,  petrochemical, 
electronics,  and  new-materials  industries  in  the  leading  position,  to  form 
a  special  characteristic  of  "having  it  when  others  lack  it,  and  having  it 
in  superior  form  when  others  have  it."  Some  suggested  using  Xiamen's 
advantages  in  oceanographic  research  to  organize  the  10-odd  research  units 
and  schools  together  to  begin  comprehensive  oceanographic  development. 

The  provincial  leaders  Xiang  Nan  [7309  0589],  Jia  Qinglin  [6328  1987  2651], 
and  Cai  Ninglin  [5591  1337  2651]  came  to  the  meeting  to  hear  the  suggestions 
given  by  everyone  for  "locomotive"  industries  in  our  province's  S&T  develop¬ 
ment  plan. 

This  project  was  begun  last  year  under  the  leadership  of  the  leaders' 
task  force  for  provincial  S&T  planning.  They  divided  into  24  specialized 
groups  and  organized  the  participation  of  more  than  1,100  famous  experts 
and  scholars  in  discussions  and  working  out  plans.  They  forecast  that  at 
the  end  of  this  year  they  will  have  worked  out  a  medium-  and  long-term  S&T 
development  plan  for  our  province. 
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MC68000  MONITOR  PROGRAM  DISASSEMBLED 

Chongqing  WEIXING  JISUANJI  [MICROCOMPUTERS]  in  Chinese  No  2,  14  Mar  85 
pp  102-107 

[Article  by  Yuan  Hongchun  [5913  1347  2504],  Chengdu  Academy  of  Electronic 
Communications  Engineering:  "A  Monitor  for  the  MC68000  Single  Board  Micro¬ 
computer"] 

[Text]  I.  Preface 

The  Chengdu  Academy  of  Electronic  Communications  Engineering  [CAECE]  Micro¬ 
computer  Development  and  Training  Center  and  the  state  run  Factory  830,  on 
the  basis  of  selling  the  foreign  MC68000  DEC  (MC68000  Educational  Computer 
Board)  sample  materials,  have  cooperatively  developed  the  successful 
DJS-0600N  16-bit  single  board  microcomputer,  which  underwent  technical 
evaluation  in  early  July  1984.  It  is  now  being  produced  in  small  quanti¬ 
ties  at  Factory  830. 

In  the  process  of  completing  this  research  work,  the  author  has  disassembled 
and  printed  out  the  16K  byte  monitor  that  is  burned  into  the  MC68000  DEC 
EPROM  on  this  single  board  computer.  Based  on  the  program  inventory  so 
obtained,  he  has  analyzed  the  operating  processes  of  the  26  monitor  commands 
in  this  monitor. 

With  the  DJS-0600N  single  board  computer  as  background,  this  paper  briefly 
discusses  the  important  contents  of  the  monitor  routines. 

The  author  feels  that  if  we  are  to  further  study  and  analyze  microcomputer 
monitor  programs ,  or  if  we  hope  to  utilize  or  develop  in  our  own  work  the 
same  ideas  in  microcomputer  software,  then  this  paper  can  provide  a  certain 
amount  of  help  and  reference  value. 

My  classmates  Su  Guoquan  and  Cheng  Dongyang  of  the  CAECE  computer  department 
joined  in  this  monitor  analysis,  doing  a  great  deal  of  work.  I  extend  my 
heartfelt  appreciation. 

Due  to  my  limitations,  there  will  be  certain  errors  and  deficiencies,  for 
which  I  ask  instruction  from  my  readers . 
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II.  Overview  of  the  DJS-0600N 

2.1  Important  Technical  Characteristics 

The  DJS-0600N  single  board  microcomputer  has  at  its  core  the  MC68000  micro¬ 
processor  chip,  which  with  other  LSI  devices  is  packaged  on  a  30  cm  x  24  cm 
two-sided  printed  circuit  board. 

Chief  technical  characteristics: 

_ 4  mHz  MC68000  16-bit  MPU  can  provide  to  the  user  eight  32-bit  data  and 

address  registers,  respectively,  a  range  of  16  megabytes  directly  address¬ 
able,  56  powerful  instructions ,  5  kinds  of  operations,  and  6  basic  types  and 
14  flexible  addressing  modes. 

— ‘Sixteen  16K  bit  dynamic  RAM's  (MCM4116),  composing  32K  bytes  of  user  memory. 

— Two  8K  byte  EPROM's  (MCM68A7 64)  making  up  16K  bytes  of  EPROM.  The  monitor 
program  called  "TUTOR"  is  burned  into  the  EPROM,  which  can  provide  monitor, 
debugging,  and  assembly /disassembly  functions. 

— 2  RS232C  compatible  serial  input /output  ports  made  up  of  MC6850ACIA's. 

—One  16-bit  parallel  input /output  interface  with  a  24-bit  programmable  timer, 
composed  of  an  MC68230  (PI/T). 

— 1  audio-frequency  tape  input/output  port. 

_ Fitted  with  RESET  and  ABORT  function  buttons ,  making  it  easier  for  the  user 

to  end  a  program  under  abnormal  conditions,  and  to  return  again  to  the  monitor 

status. 

In  use,  the  user  need  only  connect  an  external  CRT  and  a  DC  power  source 
(+5  V,  ±12  V)  to  form  a  basic  system  with  reasonable  capabilities.  It  is  a 
powerful  method  and  tool  for  the  study  and  development  of  MC68000  applications. 

2.2  Schematics 

This  schematic  shows  the  main  devices  used  in  this  single  board  computer  and 
the  logical  connections  between  them. 

We  point  out  also  at  this  time  that  on  this  single  board  computer,  port  No  1 
is  provided  for  CRT  use.  Port  No  2  is  for  modem  or  mainfrane  use.  Port  No  3 
is  for  a  line  printer.  Port  No  4  is  then  for  connecting  a  recorder  to  act  as 
a  peripheral  storage  device. 
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2.3  Storage  Area  Address  Allocation 


Function 

Address 

[Store  abnormal  vector  table 

Jrom/eprom 

$000000-$000007 

System  storage  i 

[RAM 

$000008~$0003FF 

Monitor  program  uses  scratch 
„  area 

RAM 

$000400-$0008FF 

User  storage 

RAM 

$000900-$007FFF 

Monitor  program  (TUTOR) 

Firmware  package 

ROM/EPROM 

$008000-$00BFFF 

Unused  area 

1 

[kl/T  parallel  interface /timer 

$00C000-$00FFFF 

$010000-$01003F 

I/O  equipment  < 

ACIA2  (low  byte) 

$010040-$010043 

1 

and  ACIAi  (high  byte) 

$01FFFF 

Unused  area 

$020000-$02FFFF 

MC6800  page 

$030000-$03FFFF 

Unused  area 

$040000-$FFFFFF 
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III.  Monitor  Program 

This  monitor  program  is  called  TUTOR*  is  written  in  MC68000  mnemonics,  and 
occupies  a  total  of  16K  bytes. 

TUTOR  is  chiefly  composed  of: 

1.  routines  for  system  initialization 

2.  routines  to  enter  the  monitor 

3.  main  routines  to  accomplish  26  command  functions 

4.  various  subroutines  used  by  the  main  command  routines 

5.  more  than  20  basic  common  subroutines 

6.  byte  string,  byte  table,  data,  working  units,  etc. 

7.  assembly  routines 

8.  disassembly  routines 

Implementation  of  all  routines  uses  the  structured  method  of  calling  sub¬ 
routine  during  the  process  of  running  the  main  routines ,  while  subroutines 
themselves  can  be  nested  for  calling. 

The  relations  between  TUTOR  and  the  user  are  realized  by  the  user  typing 
certain  commands  on  the  CRT  terminal  keyboard.  Those  commands  may  be  divided 
into  four  types: 

1.  allows  user  to  display  or  modify  memory 

2.  allows  user  to  display  or  modify  each  internal  register  of  the  MC68000  MPU 

3.  allows  user  to  run  programs  at  various  levels  of  control 

4.  control  access  to  all  I/O  on  the  single  board  computer 

There  is  yet  another  function  in  TUTOR,  which  is  the  so-called  "TRAP14"  han¬ 
dling  routine.  This  allows  the  user's  program  to  call  certain  subroutines  in 
TUTOR  directly. 


The  operation  mode  of  TUTOR  is  as  shown  in  the  following  diagram 


9 


On  power  up  or  by  pressing  the  RESET  button,  the  single  board  computer  will 
begin  system  initialization.  The  CPU  is  reset  and  the  TUTOR  program  is  called. 
After  initialization,  the  terminal  screen  will  print  the  following  prompt: 

TUTOR  l.X> 

(where  X  is  the  software  version  number) 

It  then  waits  for  a  response.  At  this  time,  the  user  can  then  type  in  any 
command  supported  by  the  firmware  package.  While  the  user  is  still  inputting 
the  first  line,  the  standard  input  subroutines  in  TUTOR  will  control  the 
system,  and  only  when  a  carriage  return  has  been  input,  indicating  completion 
of  the  line,  will  command  processing  begin. 

3.1  TUTOR  Initialization  Routines  and  Entering  the  Monitor 

When  the  system  has  been  initialized,  the  MC68000  microprocessor  will  auto¬ 
matically  fetch  the  beginning  addresses  of  the  system  supervisor  mode  stack 
pointer  and  initialization  routines  from  the  eight  cells  from  $8000  to  $8007, 
and  place  them  in  SS  and  PC,  respectively.  It  then  begins  execution  of  the 
initialization  routines. 

The  routine  flowchart  is  as  shown  below.  Blocks  that  have  addresses  to  the 
right  indicate  subroutines  pertaining  to  that  block,  which  address  is  then 
the  entering  address  for  the  subroutine  (the  same  with  the  rest). 


_ ,L  $8268 _ 

output  the  contents  of  the  buffer  to  .  ' 

terminal  connected  to  port  1 _  ?9C14 

_ 1  $826C _ 

initialize  buffer  pointers  A5,  A6  |  $9CA0 

_ _ 1  $8270 _ 

input  a  byte  from  the  keyboard  $9CA8 


Figure  2. 
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Command  name 

Meaning 

Entry  point 

BF 

memory  block  fetch 

$84E8 

BM 

block  move 

$8538 

BR 

breakpoint  set 

$8654 

BS 

block  search 

$871A 

BT 

block  test 

$8854 

DC 

data  conversion 

$88A6 

DF 

register  group  for  display  formatting 

$8904 

DU 

dump  to  memory 

$8A22 

GD 

direction  execution 

$8634 

GO 

run  a  program 

$8626 

GT 

run  to  a  breakpoint 

$85DA 

HE 

help  (display  all  command  names) 

$8E84 

LO 

load 

$8FC8 

MD 

memory  display 

$9166 

MM 

modify  memory 

$931E 

MS 

memory  set 

$94A2 

NOBR 

clear  breakpoints 

$9650 

NOPA 

printer  disable 

$A1C2 

OF 

display  offset  value 

$96C6 

PA 

printer  connect 

$A1C6 

PF 

port  mode 

$9744 

•  RX 

one  register  display/edit 

$9596 

TM 

transparent  mode 

$987A 

TR 

tracing 

$858A 

TT 

temporary  breakpoint  tracing 

$85A6 

VE 

verify 

$8FC2 

* 

send  data  to  port  No  2 

.A  -.A 

display/set  address  registers 

.D  -.D 

display /set  data  registers 

.R  -.R 

display/set  corresponding  shift  registers 

.PC 

display /set  program  data 

.SR 

display/set  status  registers 

.SS 

display /set  sup.  stack  pointer 

.US 

display/set  user  stack  pointer 

11 
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A 

,$82D8 

Set  system  in 

supervisor  mode 

,$82DC 

_ T  - - - - - 

Add  2  to  buffer  pointer 


1 - - - 

N 

i$82DE-$82E6 

Send  the  true  address  found  out  for  the  command  byte 
t-n  A 7.  then  send  to  Ao,  while  (A7>+4  +A7 

U - =- - - - 

V 

-S82E8 

Jump  to  the 
and  execute 
the  command 

location  for  which  (Ao)  has  the  address 
,  i.e.,  transfer  to  the  entry  point  of 
selected 

»  - - -  1  '  - - 

,$8232 

(  return  to  TUTOR 


Figure  4 

T3 


In  these  routines  the  user  input  is  checked  to  see  what  command  it  is.  which 
depends  upon  comparing  the  command  input  by  the  user  against  a  command  table 
set  up  in  a  section  of  the  program  memory.  In  this  section  of  memory 
addresses,  each  four  bytes  is  used  by  a  command.  The  first  two  bytes  store 
the  command  name  and  the  latter  two  bytes  store  the  difference  between  the 
command  entry  point  and  the  $8000  location.  When  the  user  issued  command 
corresponds  to  one  in  the  command  table,  the  routine  jumps  to  the  command 
entry  point  and  executes  the  command. 

3.2  The  Names  of  Monitor  Commands  and  Entry  Points 

The  names  and  entry  points  of  the  TUTOR  commands  are  shown  in  the  table. 

The  last  eight  commands  in  the  table,  having  to  do  with  displaying  or  setting 
registers,  are  not  just  commands  provided  with  entry  points,  but  are  derived 
from  the  .RX  command. 

3.3  Analysis  and  Examples  of  Monitor  Command  Routines 

We  already  know  that  the  principle  for  writing  any  command  routine  lies  in 
accomplishing  the  function  belonging  to  the  command,  step  by  step.  The  under¬ 
lying  design  philosophy  for  this  monitor  program  is  the  same.  In  the  TUTOR 
monitor  program,  the  command  routine  structure  uses  levels  of  subroutine 
nesting,  as  well  as  the  method  of  having  the  main  program  call  subroutines. 

Routine  flowcharts  for  two  commands  have  been  appended  to  this  article  as 
examples. 

3.3.1  The  BM  (Block  Move)  Command 

The  function  of  this  command  is  to  move  the  contents  of  a  clock  of  memory  to 
another  block. 


The  instruction  syntax  is  as  follows: 

BM  <address  1>  <address  2>  <address  3> 

where  <address  1>  is  the  first  address  for  the  source  block,  <address  2>  is 
the  last  address  of  the  source  block,  and  <address  3>  is  the  first  address 
of  the  object  block. 

3.3.2  The  TR  (Trace)  Command 

This  command  is  used  to  trace  instructions  during  execution. 

The  command  syntax  is  as  follows: 

TR  (<counter  value>) 

The  counter  value  in  this. command  is  the  number  of  instructions  to  trace  at  one 
time.  If  no  value  is  provided  with  the  command,  this  means  that  each  command 
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(  The  command  "BM") 


_ ‘  * 

?853C 

See  whether  buffer  is  empty;  if 
empty,  then  jump  to  output  of 
"SYNTAX  ERROR"  and  return  to 

TUTOR 

$8540 

Fetch  source  block  starting 
address  from  buffer,  convert  to 
hexadecimal  and  store  at  Do 

$8544 

Send  first  address  of  source 
block  to  A3 

$8546 

Send  whether  buffer  is  empty; 
if  so,  then  jump  to  output 
"SYNTAX  ERROR"  and  return  to  § 

TUTOR 

S854A 

Fetch  the  last  address  of  source 
block,  convert  to  hexadecimal  < 

and  store  at  Do 

$854E 

Send  last  address  of  source 
block  to  Ai 

$8550 

See  if  buffer  is  empty,  if  so, 
then  jump  to  output  "SYNTAX 
ERROR"  and  return  to  TUTOR 

$8554 

Fetch  first  address  of  object 
block  from  buffer,  convert  to 
hexadecimal  and  store  at  Do 

,$8B58 

Send  first  address  of  object 
block  to  A2 

$855A-$855C 

Output  "PHYSICAL  ADDRESS=" 
source  first  and  last  addresses 

$8B3C 


i$9A7E 


$8560-$8562 


Output  "PHYSICAL  ADDRESS=" 
object  block  first  address 


Is 

-^tddre!ltof 
source  less/ 

thaad^reil? 


$856A 


Length  of 
source  block 
to  Di 


Output  "SYNTAX 
ERROR"  and  re¬ 
turn  to  TUTOR 


$80B6 


$8B3C 


$9A7E 


$8B3C 


/-/first'\ 

N  /'address  of 

- -<block<first 

-address  of 
'  source 

Y 

_ _  I  $8374 

Transfer  the  contents  of 
the  first  address  of 
source  block  to  first 
address  of  object  block; 
increment  source  block 
address  and  object  block 
address _ 


/block  £.  \  N 
Contents  fin-^111 

/ished  trans- 
-  .  ?  .. 


T 


.1  $8580 


Transfer  contents  of 
source  block  last 
address  to  object  block 
last  address;  decrement 
source  and  object  block 
final  addresses _ _ 


I 

contents 
finished 
\transrer 
? 


Figure  5 


(  return  to  TUTOR  ) 
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(  The  command  "TR"  ) 

— y 

_i_  $858E 


Check  whether  buffer  is  empty; 
if  so,  set  up  as  step  tracing 


$8592 


Convert  quantity  to  be  traced 
to  hexadecimal  and  store  at  Dp 


j  $8596 


See  whether  quantity  to  be 
traced  exceeds  address  limits 


$9A7E 


$9B00 


$9AA2 


Store  1  at  Df 


L 


Quan- 

Y  tity  to 
'“*6\be  traced 
'\equal  to^ 

"  .N 


Send  quantity  to  be  traced 
(Dn)  to  $04DA 


#  -  1  to  $04DE  sets  system  in 
trace  mode 

- - - 

$84CE-$85D2 

Output  "PHYSICAL  ADDRESS="  the 
first  address  of  the  instruc¬ 
tion  to  be  traced 

$80BC 

■ 

. 

' 

Set  up  NMI  interrupt  handling 
routine  entry  point  at  $8CBC 
to  go  to  $0024 

------  -  ^ 

$8D5E 

Temporarily  store  the  monitor  program  in  the 
buffer  and  send  the  contents  of  the  stored 
registers  to  the  registers  corresponding  to 
the  MPU,  then  put  the  contents  of  PC,  SS  on 
the  stack 

'© 


Take  the  contents  of  PC*  SS  from  the  stack  and 
send  to  the  corresponding  MPU  registers,  and 
begin  execution  of  the  routine  from  the  current| 
line  number  in  the  PC 


Figure  6 
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is  to  trace  one  execution.  The  initial  address  of  the  instruction  to  be 
executed  is  given  by  the  value  of  the  program  counter  (PC) .  After  each 
instruction  has  run,  this  command  will  also  display  on  the  CRT  screen  the 
contents  of  the  microprocessor’s  internal  registers  at  this  time  (PC,  SR, 

US,  SS,  D0-D7,  A0-A7)  for  the  user  to  analyze  and  use.  At  the  same  time! 
it  will  disassemble  and  display  the  instruction  at  the  address  pointed  to 
by  the  PC. 

3.4  Examples  of  Common  Subroutines 

In  the  TUTOR  monitor  program  there  are  also  more  than  20  commonly  used  sub¬ 
routines,  as  for  example:  an  input  subroutine  that  inputs  a  byte  from  any 
port  to  the  single. board  computer,  an  output  subroutine  that  outputs  a  byte 
to  any  port,  a  conversion  subroutine  that  converts  1,  2,  4,  6,  or  8  hexadeci¬ 
mal  data  to  ASCII  code  or  that  converts  ASCII  code  to  hexadecimal  data,  as 
well  as  several  subroutines  that  accomplish  other  functions. 

In  addition  to  providing  command  routine  calls,  in  this  program  common  sub¬ 
routines  may  be  called  directly  by  the  user  by  calling  the  "TRAP14"  handling 
routine.  6 

For  an  example,  we  will  here  use  only  the  input  routine  that  sends  data  from 
the  keyboard  to  a  memory  buffer  to  analyze  the  operating  process  of  this 
subroutine. 

As  we  said  before,  we  know  that  there  are  four  ports  in  this  single  board 
computer.  Two  serial  ports  connect  a.  CRT  and  a  host  computer,  respectively, 
(or  a  modem) ,  while  the  other  two  ports  are  used  to  connect  a  line  printer 
and  a  tape  cartridge  recorder.  Therefore,  there  are  four  I/O  subroutines. 

The  entry  points  to  these  subroutines  are  as  in  the  following  table: 


Port  No 

Routine  entry  point  stored  at 

Routine  entry  point 

I/O 

Port  1 

$062C 

$9CC8 

input 

$0630 

$9C14 

output 

Port  2 

$0634 

$9FDE 

input 

$0638 

$9C34 

output 

Port  3 

$06  3C 

$9CC8 

input 

$0640 

$9D70 

output 

Port  4 

$0644 

$A0A6 

input 

$0648 

$9EA6 

output 

12586/9365 
CSO:  4008/1073 
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[Text]  In  this  paper,  microprocessor  software  is  used  to  correct  error  codes. 
Hamming  code  and  cyclic  code  are  interwoven  to  correct  8  bits  of  burst  error. 
Thus ,  the  use  of  complicated  equipment  required  in  the  conventional  inter¬ 
weaving  method  can  be  avoided. 


Hamming  Code 


A  universal  monitoring  matrix  H  is  chosen  to  be  a  (7,4)  Hamming  code 


11110  0 
0  110  10 
10  10  0  1 


The  code  length  is  7.  There  are  4  information  bits  and  3  monitoring  bits. 
It  is  capable  of  correcting  any  single  random  error  in  a  string  of  codes. 
Let  us  assume  that  the  information  sequence  is  C0,  Ci ,  C2,  C3  and  the  moni¬ 
toring  sequence  is  C* ,  C5,  C6 ,  then  the  code  block  (C0 ,  Ci,  C2 ,  C3,  C4 ,  C5, 
Cg)  satisfies  the  following: 


r  0  1 

1  o 
L 1  1 


i.e.,  HCT  =  0T 

where  0  =  [0  0  0] 


1110 
1  l  0  1 

0  10  0 


On 

0 

1 


/ 

C, 

c, 

Cs 

c 4 

C6 


r0l 

0 

10  J 


(1) 
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It  can  be  expressed  by  the  following  equations: 


(  C{~Cl§)Ci®Cs 

jc3=c0©c2©c3 

(ce  =  6’o©C,©C3 


(3) 


Let  us  assume  that  the  error  pattern  created  when  the  code  block  is  transmitted 
in  a  channel  is  E  =  (E0E1E2E3E4E5E6).  If  there  is  an  error  in  the  ith  bit, 
then  E^  =  1,  otherwise,  it  is  zero.  If  the  code  block  received  is 
C'  =  C+E  =  (Co  C{  C2  C3  C4  C5  Ce),  where  +  E^,  based  on  equation  (2) 

we  get 


Let  ST 
(5) 


HC,t  =  H  (C+E)T  =  HCT  +  HEt  =  HEt  (4) 

HeT,  or  S  =  EH^,  where  S  is  the  conjugate  transpose.  We  get  equation 


■  5,' 

*0 

1 

1 

1 

1 

0 

o' 

= 

1 

0 

1 

1 

0 

1 

0 

.  S3_ 

.1 

1 

0 

1 

0 

0 

1 . 

(Et\ 

E> 

E1 

Et 

Et 

E, 


(5) 


In  transmission,  if  there  is  only  an  error  at  C2 ,  then  E  =  (0010000).  Substi¬ 
tuting  it  into  equation  (5),  we  get  S  =  [S2SiSol  =  [110].  Therefore,  we  can 
see  that  if  only  Ej  is  1  in  E,  ST  is  the  ith  column  in  the  matrix  H.  Table  1 
is  tabulated  based  on  equation  (5). 


In  equation  (4),  ST  =  HC,T,  i.e. 


'  Si  1 

'0 

] 

1 

1 

1 

0 

o’ 

= 

1 

0 

1 

1 

0 

1 

0 

.  S3_ 

.  1 

1 

0 

1 

0 

0 

1 . 

jC'0\ 

ci 

Ct 

cl 

Cl 

Cl 

[Cl) 


Then,  the  decoding  equations  are: 

fS2  =  C'©C£©C£©C4' 

j^^©^©^©^  (7 

U3-C6©C{©C'S©CS 

By  using  the  coding  equations  in  (3),  decoding  equations  in  (7)  and  Table  1, 
error  correction  can  be  carried  out  by  software. 
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Table  1.  Relation  Between  Error  Location  and  Conjugate  Transpose 


Error  location 

S 

Error  location 

S 

S2 

Si 

So 

S2 

Si 

So 

0 

0 

0 

c3 

1 

1 

0 

Co 

0 

1 

1 

Ci* 

1 

0 

0 

Ci 

mm 

0 

1 

C5 

0 

1 

0 

C2 

n 

1 

0 

C6 

0 

0 

1 

Because  there  are  numerous  I/O  interfaces  and  transmission  modes,  we  will  not 
discuss  this  issue  in  this  paper.  We  will  only  introduce  code  compiling  pro™ 
cess  using  a  Z80  microprocessor  or  a  single  board  computer  software.  In 
order  to  accelerate  the  code  compiling  rate,  assembly  language  is  used. 

1.  Encoding 

In  Figure  1,  the  front  half  byte  (a0-a3)  of  the  first  address  of  the  input 
data  is  coded  and  transmitted  to  the  first  output  address.  Then  the  latter 
half  byte  (a4-a7)  is  coded  and  sent  to  the  second  byte.  The  same  process  is 
repeated  for  the  second  input  byte  until  all  the  data  are  coded.  The  coding 
process  is  shown  in  Figure  2.  In  equation  (3),  because  the  value  of  Ci  in  the 
information  and  monitoring  bits  is  in  the  GF92)  domain,  the  addition  of  the 
module  2  of  Ci,  C2,  and  C3  is  equivalent  to  the  ordinary  addition  among  them¬ 
selves;  which  is  1  for  odd  numbers  and  0  for  even  numbers.  The  latter  is 
adopted  to  simplify  the  program  significantly.  In  the  program,  registers  C, 

D,  and  E  are  assigned  to  represent  the  values  of  the  monitoring  bits  C4  ,  C5, 
and  C6 ,  respectively. 


Input 


Figure  1. 


In  Figure  2,  RRA  is  the  cyclic  right  shift  of  the  contents  of  the  accumulator 
A  and  the  carry  marker.  From  the  carry  marker,  we  can  determine  whether  the 
Oth  bit  of  A,  Co,  is  0  of  1.  If  it  is  1,  then  B0  is  set  to  be  1.  Otherwise, 
Bo  is  set  to  be  0.  B0  is  the  Oth  bit  in  the  B  register.  Based  on  equation 
(3),  if  C0  is  1,  then  1  is  added  to  both  registers  D  and  E.  The  immediate  RRA 
to  follow  is  to  process  Cx.  After  four  such  RRA,  the  content  in  register  A, 

C0  -  C3,  is  first  entered  into  a0  -  a3  in  register  B  and  then  used  to  evaluate 
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registers  C,  D,  and  E  in  order  to  ascertain  the  contents  in  monitoring  bits 
Ci*  >  C5 ,  and  Cg  and  send  them  to  a4-a7  in  B.  The  next  cycle  is  to  code  the 
latter  half  byte  of  the  input  data  and  send  it  to  the  second  byte  of  the  out¬ 
put  data  until  the  entire  coding  process  is  completed.  In  Figure  2,  the 
address  of  the  first  data  block  should  be  an  even  number  to  allow  the  program 
to  proceed  as  shown  in  Figure  1. 
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2.  Decoding 

Based  on  equation  (7),  decoding  is  done  at  the  receiving  end.  We  only  correct 
errors  in  information  bits  Co  -  C3  and  ignore  monitoring  bit  errors.  The 
decoding  process  is  schematically  shown  in  Figure  3.  Table  2  is  the  Hamming 
code  error  correction  table.  If  an  error  occurs  at  C2 »  its  conjugate  S  =  (110) 
appears  as  06H  in  a  register.  This  number  is  called  the  conjugate  number. 

Its  corresponding  error  pattern  is  E2  =  1  and  the  remaining  =  0.  In  the 

corresponding  register,  an  error  occurs  at  C2  and  the  register  value  is  04H, 
which  is  called  the  correction  number.  In  the  program,  04H  is  used  to  correct 
the  error  in  C2.  In  Figure  4  the  symbol  CORREC  is  used  to  represent  the  first 
addresses  of  the  conjugate  number  and  the  correction  number  (CORREC:  DB03,  05, 
06,  07,  01,  02,  03,  04,  08). 


Input 


Figure  3. 


Table  2. 


Error  location 

Conj ugate 
S2  Si 

So 

Conjugate  number 

Error  correction  number 

f 

Co 

0 

1 

1 

03H 

01H 

I 

Ci 

1 

0 

1 

05H 

02H 

1 

c2 

1 

1 

0 

06H 

04H 

1 

c3 

1 

1 

1 

07H  7.. 

08H 

In  Figure  4  the  values  of  S2,  Si,  and  So  in  equation  (7)  are  stored  in  regis¬ 
ters  B,  C,  and  D,  respectively,  from  the  start  to  the  point  W  where  decoding 
of  equation  (1)  is  completed.  The  conjugate  transpose  is  obtained  atfW.  It 
is  compared  to  conjugate  number  of  the  first  address  of  CORREC.  If  C2  is 
wrong,  the  conjugate  transpose  in  A  is  110,  which  is  equal  to  the  conjugate 
number  06  in  HL.  Then,  HL  plus  3  is  04  (the  correction  number).  The  error 
in  C2  is  corrected  by  using  the  XORB  command.  The  following  program  is  used 
to  complete  the  program  shown  in  Figure  3  with  the  first  address  of  the  data 
block  being  even.  This  process  continues  until  data  processing  is  completed. 

Cyclic  Code 

Let  us  choose  the  generating  polynomial  g(x)  to  be: 

g(x)  =  x3  +  x2  +  1  (®) 
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Figure  4. 
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Then,  the  generating  matrix  G  is: 


-1000110 
0  1  0  0  0  1  1 

0  0  1  0  1  1  I 

Lo  0  0  1  1  0  1J 


(9) 


Thus,  a  (7,4)  cyclic  code  for  correcting  a  1  bit  random  error  can  be  totally 
defined.  In  order  to  facilitate  encoding  and  decoding,  it  is  necessary  to 
make  the  generating  matrix  G  be  a  universal  monitoring  matrix  H: 


0 

1 


Lo  l 


11100' 
1  0  0  1  0 
11001. 


(10) 


The  encoding  equations  are : 

c^c.ec.ec, 

c5=c„®c,®cs  (11) 

iC<  =  C1®C;®C3 


and  the  decoding  equations  are 

S,  =  C(>®C{®C£®Cj  (12) 

.5',-C{®C5®C3,®C,/ 

Equation  (11)  is  similar  to  (3)  and  equation  (12)  is  similar  to  (7),  therefore, 
the  code  compiling  program  is  also  similar. 


1.  Encoding 

The  encoding  scheme  is  shown  in  Figure  1.  Figure  5  shows  the  cyclic  code 
compiling  program.  Registers  C,  D,  and  E  are  used  to  represent  the  values  of 
C4 ,  C5 ,  and  Cg  in  equation  (11),  respectively. 

2.  Decoding 

The  decoding  scheme  is  shown  in  Figure  3.  Figure  6  shows  the  flow  chart  of 
the  cyclic  code  decoding  program.  Registers  B,  C,  and  E  represent  S2,  Si,  and 
S0  in  equation  (12),  respectively.  The  error  correction  table  is  shown  in 
Table  3,  using  CORREC  to  represent  the  conjugate  number  and  the  first  correc¬ 
tion  number  address  (CORREC:  DB06,  03,  07,  05,  01,  02,  04,  08). 


Table  3.  Cyclic  Code  Correction  Table 


Error  location 

Conjugate 
S2  Si 

S0 

Con j  ugate  number 

Error  correction  number 

f 

Co 

1 

1 

0 

06H 

01H 

cl 

0 

1 

1 

03H 

02H 

cl 

1 

1 

1 

07H 

04H 

cl 

1 

0 

1 

05H 

08H 
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Interwoven  Hamming  Code 

As  shown  in  Figure  7,  at  the  sending  end,  the  first  half  input  data  byte 
Cio  -  Ci 3  is  Hamming  coded  to  Cio  -  Cig  and  transmitted  to  the  9th  bit  of 
each  corresponding  byte  in  the  first  output  data  block  (each  data  block  con¬ 
tains  7  bytes).  Then,  the  latter  half  of  the  first  input  byte,  C20  -  C23,  is 
similarly  treated  until  4  input  bytes  are  processed  and  sent  to  fill  the  first 
output  data  block.  Then,  the  next  block  is  similarly  processed.  This  is  the 
encoding  process  of  the  interwoven  code  at  the  sending  end. 

At  the  receiving  end,  the  Oth  bit  of  each  byte  in  the  first  input  block  is 
decoded  into  Ci'o  -  CiV .  After  correcting  the  errors,  Ci'o  -  Ci'3  is  sent  to 
the  front  half  of  the  first  byte,  as  shown  in  Figure  8.  It  proceeds  in  this 
manner  until  the  first  input  block  is  Hamming  decoded  and  sent  to  the  output. 
Then  the  next  block  is  processed.  This  is  the  interwoven  decoding  process  at 
the  receiving  end. 

Based  on  the  above,  8  (7,4)  Hamming  codes  form  an  interwoven  code  matrix.  The 
degree  of  interweaving  is  8.  Each  byte  contains  the  first  bit  of  every  (7,4) 
code.  In  the  transmission  process,  if  the  burst  length  is  8  bits,  then  it  is 
equivalent  to  having  a  1  bit  error  in  each  (7,4)  code  at  the  receiving  end 
during  decoding.  Each  (7,4)  code,  however,  can  correct  1  bit  of  error.  Thus, 
these  8  bits  of  burst  errors  can  also  be  corrected. 

1.  Interwoven  Encoding 

Figure  9  is  the  encoding  flow  chart.  Equation  (3)  is  used  to  encode.  Let 
registers  C,  D,  and  E  represent  the  values  of  Ci* ,  C5,  and  Ce ,  respectively. 
Shifting  (HL)  to  the  right  by  1  bit  is  done  by  using  the  RRC(HL)  command. 

Based  on  the  scheme  shown  in  Figure  7,  a  block  is  encoded  before  the  next 
block  is  processed  until  the  entire  job  is  completed. 

2.  Interwoven  Decoding 

Figure  10  is  the  decoding  flow  chart.  Equation  (7)  is  used  for  decoding.  Let 
registers  B,  C,  and  D  represent  the  values  of  S2 ,  Si,  and  So,  respectively. 
With  the  exception  that  its  error  correction  numbers  are  10,  20,  30,  and  40, 
the  rest  of  the  error  correction  table  is  the  same  as  Table  2.  Therefore, 
CORREC:  DB03,  05,  06,  07,  10,  20,  40,  and  80. 

The  Cyclic  Interwoven  Code 

The  interweaving  process  of  the  cyclic  interwoven  code  is  the  same  as  those 
in  Figures  7  and  8.  Encoding  and  decoding,  however,  are  done  based  on  equa¬ 
tions  (11)  and  (12).  Figures  11  and  12  are  the  flow  charts  for  encoding  and 
decoding  cyclic  interwoven  codes.  The  registers  are  set  as  that  for  the  Ham¬ 
ming  interwoven  code.  With  the  exception  that  the  error  correction  numbers 
are  10,  20,  40,  and  80,  the  rest  of  the  error  correction  table  is  the  same  as 
Table  3.  Therefore,  CORREC:  DB06 ,  03,  07,  05,  10,  20,  40,  and  80. 
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Based  on  the  above,  the  main  portion  between  the  cyclic  code  and  the  Hamming 
code,  as  well  as  between  the  cyclic  interwoven  code  and  Hamming  interwoven 
code,  is  identical.  Therefore,  it  is  possible  to  set  up  subroutines  to  handle 
them. 

Conclusions 

In  the  coding  process,  we  should  first  consider  fast  speed,  i.e.  less  T  states. 
Then,  we  should  conserve  the  internal  storage  elements.  For  example,  the 
command  XORA  is  three  state  number  less  than  that  of  LD  A,  00H.  LDBC,  nhH  is 
4  T  less  than  that  of  LDB ,  nH  and  LD  C,  nH.  RRA  is  4  T  less  than  that  of 

BITO,  A.  JP  is  2  T  less  than  JR,  however,  there  is  one  extra  internal  storage 

element. 

The  data  processing  rates  of  common  data  transmitters  and  demodulators  are 

I. 2  kb/s,  2.4  kb/s  and  4.8  kb/s.  Using  4.8  kb/s  as  an  example,  it  takes 

II. 667  ms  to  transmit  a  block  of  data  (a  total  of  7  bytes  in  56  bits).  The 

execution  time  required  for  a  2  MHz  microprocessor  using  a  encoding  program 
with  4,522  states  is  2.261  ms.  The  execution  time  of  a  decoding  program  with 
6,325  states  is  3.163  ms.  If  the  primary  frequency  of  the  microprocessor  is 
4  MHz,  then  the  execution  time  could  be  cut  in  half.  Thus,  at  a  data  trans¬ 
mission  rate  of  4.8  kb/s,  microprocessor  coding  has  no  effect.  It  can  be 
considered  transparent .  An  interrupt  mode  may  be  used  at  both  the  receiving 
and  sending  end  to  perform  data  input  and  output  simultaneously. 
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[Text]  Abstract:  Because  of  its  high  temporal  and 
spatial  resolution,  wide  spectral  and  dynamic  input 
intensity  range,  high  gain,  and  ability  to  record  tran¬ 
sient  single  events,  the  proximity  focused  microchannel 
plate  (MCP)  image  intensifier  can  be  used  as  a  super 
high-speed  optical  shutter. 

A  microchannel  plate  image  intensifier  and  a  gating 
pulse  generator,  which  controls  its  shutter  effect, 
constitute  a  super  high-speed  camera  with  high  gain 
and  fast  response .  In  this  paper ,  a  comparison  of 
three  pulse  control  modes  for  a  typical  second  genera¬ 
tion  image  intensifier  (18  mm.  Model  ITT-F4111  image 
intensifier)  is  made.  Based  on  the  comparison,  a 
novel  control  circuit  for  gating  the  photocathode  of 
the  image  intensifier,  as  well  as  capable  of  generat¬ 
ing  high  amplitude  nanosecond  pulses,  is  introduced. 

It  was  experimentally  determined  that  the  exposure 
time  of  the  camera  could  reach  approximately  2  ns  by 
using  this  control  generator.  In  addition,  the  paper 
also  discusses  the  match  between  the  generator  and  the 
image  intensifier.  Furthermore,  a  simple  and  effec¬ 
tive  method  is  introduced  to  eliminate  multiple  expo¬ 
sures  which  must  absolutely  be  avoided  in  single  shot 
photography. 

I.  Introduction  of  the  Problem 

The  proximity  focused  microchannel  plate  image  intensifier  is  a  novel  optoelec¬ 
tronic  device.  In  addition  to  being  used  as  an  optical  shutter  alone,  it  also 
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has  a  significant  light  intensifying  effect.1  The  combination  of  a  micro- 
channel  plate  image  intensifier  shutter  and  high-speed,  high-voltage  pulse 
generator  to  control  the  optical  shutter  can  constitute  a  super  high-speed 
camera;  capable  of  capturing  single  instantaneous  events  with  an  extremely 
short  exposure  time  and  high  optical  gain.  This  camera  has  a  wide  range  of 
applications ,  involving  many  important  technical  areas  such  as  high-voltage 
gas  discharge  and  laser  nuclear  fusion 

Without  any  doubt,  the  performance  of  the  camera  is  mainly  determined  by  the 
image  converter  tube  itself.  There  is,  however,  the  problem  of  how  to  more 
effectively  utilize  the  image  converter  tube.  Particularly,  two  major  per¬ 
formance  specifications — exposure  time  and  optical  gain — of  the  image 
converter  tube  of  the  microchannel  plate  image  intensifier  shutter  depend 
directly  on  the  width  and  amplitude  of  the  control  pulse.  Therefore,  the 
generation  of  large  amplitude  narrow  pulses  to  control  the  optical  shutter 
is  a  key  technical  issue  in  the  development  of  an  image  converter  tube  for  a 
super  high-speed  shutter  camera. 

In  this  paper,  a  typical  second  generation  image  intensifier,  the  Model  ITT- 
F4111  proximity  focused  microchannel  plate  image  intensifier,  is  used  as  an 
example  to  illustrate  the  pulse  control  of  a  high-speed  image  intensifier 
optical  shutter. 

II.  Three  Control  Modes  for  an  Image  Intensifier 

The  Model  ITT— F4111  image  intensifier  is  a  compact  specialized  image  converter 
tube  and  its  external  dimensions  are  shown  in  Figure  l.2  The  image  intensi¬ 
fier  is  a  circular  plate  whose  external  diameter  is  45  mm.  The  effective 
diameter  of  the  photocathode  and  the  optical  screen  is  18  mm.  The  total 
thickness  is  approximately  22  mm. 


Figure  1.  External  Dimensions  of  the  ITT-F4111  Image  Intensifier 

In  order  to  control  the  operation  of  this  image  intensifier,  the  following 
three  voltages  are  necessary: 

(1)  the  cathode  voltage  Uc  between  the  photocathode  (terminal  1)  and  the 
input  end  of  the  microchannel  plate  (terminal  2);  (2)  the  microchannel  plate 

voltage  %cp  between  the  input  of  the  microchannel  plate  and  its  output 
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(terminal  3);  and  (3)  the  anode  voltage  UA  between  the  output  of  the  micro- 
channel  plate  and  the  anode  (optical  screen,  terminal  4). 

When  the  image  intensifier  is  used  as  a  high-speed  optical  shutter,  one  of  the 
three  voltages  mentioned  above  must  be  supplied  as  pulses,  while  the  other  two 
remain  as  dc  voltages.  Thus,  we  have  three  pulse  control  modes: 

(1)  applying  control  pulses  between  the  photocathode  and  the  microchannel 
plate  input;  (2)  applying  control  pulses  between  the  input  and  the  output 
of  the  microchannel  plate;  and  (3)  applying  control  pulses  between  the  output 
of  the  microchannel  plate  and  the  anode. 

Because  of  varying  resistance  between  electrodes  and  also  because  pulses  of 
different  amplitudes  and  polarities  may  be  required,  the  difficulties 
involved  in  implementing  these  control  modes  are  also  different. 

The  input  resistance  %  and  input  characteristics  C±  between  various  electrodes 
of  the  ITT-F4111  image  intensifier  are  measured  to  be 

between  photocathode  and  microchannel  plate  input:  IL  =  4  x  1010^,  =  23pF; 

between  input  and  output  of  the  microchannel  plate:  R±  =  3.3  x  108ft»  72pF; 

between  microchannel  plate  output  and  anode:  ^  *  3.3  x  1010ft,  C±=  7pF. 

Based  on  the  data,  the  resistance  and  capacitance  between  various  electrodes 
of  the  image  intensifier  are  very  high.  Even  if  a  nearly  ideal  rectangular 
pulse  can  be  obtained  from  a  purely  resistive  low  impedance  circuit ,  it  will 
be  seriously  distorted  after  it  is  impressed  across  the  relevant  electrodes 
of  image  intensifier.  In  this  case,  not  only  the  desired  short  exposure  time 
cannot  be  obtained  because  the  edge  of  the  control  pulse  is  slowed  down  by 
the  large  time  constant  of  the  input  end  of  the  image  intensifier,  but  also, 
due  to  resonance  caused  by  the  distribution  of  parameters ,  the  optical  gate 
is  opened  many  times,  leading  to  multiple  exposures. 

Table  1  shows  a  comparison  of  the  three  pulse  control  modes  of  the  ITT-F4111 
image  intensifier. 


Table  1.  Comparison  of  Three  Control  Modes  of  the  ITT-F4111  Image  Intensifier 


Pulse  ampli¬ 
tude  to  reach  Polarity 
maximum  of  control 

Capacitance 

optical  gain  pulse 

Ground 

load  (pF) 

Usage 

Control  pulse 
photocathode 

at 

180  negative 

microchannel 
plate  input 

23 

most 

common 

Control  pulse 

at 

800  negative 

microchannel 

72 

usually 

microchannel 

plate  output 

not  used 

plate  input 

Control  pulse 
photoanode 

at 

5200  positive 

microchannel 
plate  output 

7 

sometimes 

used 
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Figure  2.  Relative  Optical  Gain  VL  vs  Photocathode  Voltage  Ur  of  the 
ITT-F4111  Image  Intensifier 

Obviously ,  the  second  control  mode  is  usually  not  adopted  due  to  the  high 
capacitance  load.  The  third  control  mode  requires  the  largest  amplitude  for 
the  control  pulse.  It  is  relatively  difficult  to  generate  and  transmit  high¬ 
speed,  high-voltage  pulses.  This  method  is  still  used  sometimes,  however, 
because  the  capacitance  load  of  this  control  mode  is  small.  In  an  ITT-F4111 
image  intensifier,  the  variation  of  the  relative  optical  gain  VL  with  cathode 
voltage  is  shown  in  Figure  2. 2  Based  on  this  figure,  several  dozen  volts 
of  voltage  change  at  the  cathode  is  sufficient  to  cause  the  optical  gain  to 
vary  by  several  magnitudes.  Only  a  180V  voltage  pulse  is  required  to  reach 
the  maximum  optical  gain  (i.e.  Vl  =  1).  In  other  words,  it  is  most  sensitive 
to  control  the  photocathode.  Because  the  first  mode  requires  pulses  of  much 
lower  amplitude  as  compared  to  those  of  the  other  two  modes,  it  is  thus  possi¬ 
ble  to  make  the  control  pulse  narrower.  Therefore,  it  may  be  possible  to  get 
the  shortest  exposure  time.  It  is  the  most  commonly  used  control  mode. 

III.  Generation  of  Large  Amplitude  Narrow  Pulses  To  Control  the  Photocathode 
of  the  Image  Intensifier 

1.  Limitation  of  the  Basic  Avalanche  Transistor  Trigger  Circuit 

In  order  to  turn  the  optical  gating  circuit  into  solid  state,  usually  an 
avalanche  transistor  is  used  to  generate  large  amplitude  pulses  with  steep 
leading  edges.  It  is  very  difficult,  however,  to  obtain  the  desired  elec¬ 
trical  control  pulses  for  the  optical  gate  based  on  the  basic  circuit  compris¬ 
ing  a  single  transistor,  as  shown  in  Figure  3,  due  to  the  following 
reasons. 


+y» 


Figure  3.  The  Basic  Avalanche  Transistor  Trigger  Circuit 
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First ,  the  wave  form  of  the  output  pulse  of  the  basic  circuit  is  strongly 
dependent  on  the  discharge  time  constant  of  the  avalanche  tube*  Ta = (Cg  +  CL) 
Rl.  We  selected  a  relatively  good  Model  2N5551  avalanche  tube  in  our  experi¬ 
ment.  All  the  data  reported  in  this  paper  were  obtained  with  this  type  of 
tube.  Figure  4  shows  the  output  pulse  wave  forms  at  various  capacitance 
values  Cs  (RL  =  50ft).  We  can  see  that  the  amplitude  of  the  pulse  increases 
with  increasing  Cg.  The  rise  and  fall  time  of  the  pulse*  however,  also 
increase.  Furthermore,  Cg  has  a  larger  effect  on  the  trailing  edge  of  the 
pulse  than  on  the  leading  edge.  When  Cg  reaches  40  pF ,  the  pulse  amplitude 
no  longer  increases  noticeably  with  Cg.  The  rise  and  fall  time,  however, 
continue  to  increase.  Therefore,  in  order  to  get  pulses  of  steeper  edges, 
it  is  preferred  to  choose  a  small  Cg.  The  price  to  pay,  however,  is  to  sacri 
fice  some  amplitude.  When  Cg  =  36  pF ,  pulses  at  70—100  V  can  be  obtained  by 
using  a  selected  2N5551  tube  on  a  50ft  load.  It  is  obvious  that  it  is  impossi 
ble  to  get  the  pulse  amplitude  required  to  control  an  optical  gate  with  a 
basic  circuit  consisting  of  one  avalanche  tube. 


Figure  4.  Output  Pulse  Wave  Forms  of  the  Basic  Circuit  at  Various  Cg's 
(10,  20,  30,  40  pF)  at  ^  =  50ft;  Y:20V/div;  X:2  ns/div 


Figure  5.  Output  Pulse  Wave  Forms  of  the  Basic  Circuit  at  Various  CL's 
(0,  40  pF)  at  Cg  =  20  pF  and  RL  =  50ft;  Y:20V/div;  X:2  ns/div 
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Next,  the  pulse  wave  form  of  the  basic  circuit  is  easily  affected  by  capaci¬ 
tance  load  because  the  load  capacitance  CL  will  directly  affect  the  discharge 
time  constant  of  the  avalanche  tube.  Figure  5  shows  the  effect  of  the  load 
capacitance  Cj  on  the  output  wave  form.  Based  on  this  figure,  when  CL  =  40  pF, 
the  pulse  amplitude  drops  to  two- thirds  of  that  at  zero  capacitive  load. 
Furthermore,  the  half  height  width  of  the  pulse  increases  to  4  ns.  It  is 
obvious  that  the  basic  circuit  shown  in  Figure  3  cannot  be  directly  used  to 
drive  an  image  intens if ier  with  a  high  input  capacitance. 

In  order  to  eliminate  these  two  disadvantages,  a  new  combination  circuit  is 
introduced.  Through  the  use  of  two  avalanche  tubes  operating  in  parallel,  the 
avalanche  currents  are  superpositioned  on  a  common  load  to  result  in  a  pulse 
almost  twice  in  amplitude  as  that  obtained  in  the  basic  circuit  (corresponding 
to  the  same  Cg) .  The  difficulty  of  choosing  an  appropriate  avalanche  tube  is 
thus  greatly  reduced.  In  order  to  minimize  the  effect  of  the  capacitance 
between  the  electrodes  of  the  image  intens if ier  on  the  pulse  generated  by  the 
control  pulse  generator,  a  special  pulse  shaping  network  consisting  of  step 
diodes  is  introduced. 

2.  Basic  Principle  of  the  Combination  Circuit 

Figure  6  shows  the  schematic  diagram  of  the  control  pulse  generator.  Through 
the  simultaneous  triggering  of  two  avalanche  tubes,  their  currents  are  super¬ 
imposed  on  a  common  equivalent  load.  Thus,  the  resulting  pulse  amplitude  is 
almost  twice  that  of  a  single  tube.  The  leading  edge  of  the  output  pulse  is 
determined  by  the  difference  between  the  rise  time  of  the  two  tubes  and  the 
switching  time.  The  trailing  edge  of  the  pulse  is  mainly  determined  by  the 
stepping  time  Tf  of  the  step  diode  used  for  pulse  shaping.  By  adjusting  the 
positive  current  If  of  the  step  diode  and  changing  the  storage  time  Tg,  we  can 
obtain  almost  symmetric  sharp  pulses.  The  basic  principle  of  the  combination 
circuit  is  explained  in  the  following  based  on  the  diagram  shown  in  Figure  7. 
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Figure  6.  Diagram  of  the  Control  Pulse  Generator 

In  the  static  state  (without  trigger  signal  input),  diodes  Di  and  D2  are 
biased  at  +Ui  and  -U2  so  that  D2  is  on  and  D2  is  off.  The  forward  current 
If  passing  through  Dx  can  be  adjusted  by. the  bias  +Ui : 

±CV 


Rdi  +  Rdi  +  Rs  +  r/ 


(1) 
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Figure  7.  Principle  of  the  Control  Pulse  Generator 


where  rf  is  the  forward  resistance  of  Di  and  Rg  is  the  serial  resistance  of 
Di.  The  circuit  on  the  right  side  of  the  dotted  line  AA'  in  Figure  7  can  be 
expressed  by  an  equivalent  resistance  Rg 


Re  -  Rdx/KRs  +  r,  +  Rod  (  ' 

Thus,  an  equivalent  circuit  of  the  generator  such  as  the  one  shown  in  Figure  8 
is  obtained. 


Figure  8.  Equivalent  Circuit  of  the  Generator 

The  energy  storage  capacitors  Ci  and  C2  in  Figure  8  are  charged  by  the  power 
source  +UB  through  a  large  resistance  to  a  voltage  slightly  lower  than  the 
breakdown  voltage  BVCB0  of  the  transistor  across  the  collector  and  the  base. 

The  prerequisite  condition  is  that  the  breakdown  voltage  of  various  transis¬ 
tors  cannot  vary  too  much.  This  can  be  ensured  by  adding  a  small  resistor 
between  the  base  and  the  collector  of  each  transistor  and  by  choosing  the 
appropriate  transistors.  The  two  avalanche  tubes  are  expressed  by  two  ideal 
switches.  Si  and  S2,  respectively. 

If  both  avalanche  tubes  are  simultaneously  triggered  by  a  positive  pulse,  then 
switches  Si  and  S2  in  Figure  8  will  close.  Capacitors  Ci  and  C2  will  discharge 
through  switches  Si  and  S2  and  the  equivalent  resistance  Rg.  A  voltage  pulse 
Vtitj  =  (ii  +  i2)  Re  g®nerated  on  Rg.  Of  course,  the  rise  time  of  the  pu  se 
generator  on  R~  is  equal  to  the  avalanche  rise  time  Ta  of  a  single  transistor 
only  when  both  tubes  have  the  same  characteristics.  If  the  two  avalanche  tubes 


39 


have  different  rise  times  or  there  is  a  difference  in  their  on  times,  the 
leading  edge  of  the  composite  pulse  will  be  slowed  down.  A  "twist"  effect, 
such  as  the  one  shown  in  Figure  9,  may  even  result.  Therefore,  we  must  select 
the  avalanche  tubes  used  in  the  combination  circuit. 


(a) 

ideal  case 


Figure  9.  Superposition  of  Two  Avalanche  Pulses 


Now,  let  us  explain  the  operating  principle  of  the  circuit  on  the  right  side 
of  the  dotted  line  AA'  in  Figure  7.  As  described  before,  in  the  static 
state,  the  diode  Di  is  biased  forward  and  its  forward  resistance  r^  is  very 
small.  Therefore,  points  ®  and  ®  in  the  figure  can  be  considered  as  a 
short  circuit.  The  diode  D2 ,  on  the  other  hand,  is  cut  off  in  the  inverse 
direction.  When  both  avalanche  tubes  are  Simultaneously  triggered,  a  large 
amplitude  negative  pulse  will  appear  near  point  @  .  Due  to  the  charge  storage 
effect  of  the  step  diode,  the  diode  Di  is  still  conducting  within  the  storage 
period  Ts  which  is  determined  by  the  forward  current  If.  Therefore,  the  nega¬ 
tive  pulse  will  also  reach  point  ®  .  At  this  time,  the  reverse  current  Ir  is 
the  step  diode  Di  will  continue  to  flow  until  the  minority  carriers  stored  near 
the  junction  is  totally  depleted.  After  the  storage  period  Ts,  the  diode  Di 
will  rapidly  enter  a  reverse  cutoff  state.  Thus,  the  reverse  current  will 
steeply  reach  its  saturation  limit  (within  the  stepping  time  Tf).  Figure  10 
shows  the  pulse  shaping  effect  of  the  step  diode  on  the  trailing  edge  of  the 
pulse.  Figure  10(a)  shows  an  avalanche  pulse  prior  to  being  shaped.  Its 
trailing  edge  is  much  slower  than  its  leading  edge.  Figure  10(b)  shows  the 
current  wave  forms  corresponding  to  diodes  with  various  Tg.  In  the  figure. 
Curve  2  is  the  optimum  case  where  the  storage  time  Tg  is  over  once  the  ava¬ 
lanche  pulse  amplitude  reaches  its  maximum.  In  this  case,  the  amplitude  of 
the  output  pulses  ,  the  maximum  and  the  base  width  of  the  pulse  is  equal  to  the 
sum  of  the  avalanche  rise  time  Ta  and  the  step  time  Tf .  The  following  rela¬ 
tion  exists  for  Ts,  If  and  Ir3 : 

Ts  -  rln(l  +  ji)  (3) 

where  T  is  the  lifetime  of  the  minority  carrier.  Because  it  involves  the 
shaping  of  large  amplitude  pulses,  the  reverse  current  Ir  in  equation  (3)  is 
very  large.  In  order  to  have  a  certain  storage  time  Tg  (which  is  necessary 
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Figure  10.  Using  Step  Tube  To  Improve  the  Trailing  Edge  of  the  Pulse 

to  shape  the  trailing  edge  of  the  pulse) ,  we  have  to  increase  If  or  T.  In 
addition,  the  diode  Di  must  have  a  high  breakdown  voltage  and  be  able  to 
sustain  high  power  pulses.  This  is  actually  too  critical  a  requirement  for 
the  step  diode.  To  this  end,  the  step  diode  used  for  pulse  shaping  has  been 
specially  designed.4  Suitable  tubes  may  also  be  chosen  from  high  reverse 
voltage  microwave  variable  capacitance  step  power  tubes. 

The  high-speed  Schottky  diode  D2  is  used  to  further  shape  the  pulses.  A 
reverse  bias  UD2  is  applied  on  the  diode  D2  to  regulate  the  adjustment  of 
— U2  in  order  to  eliminate  the  oscillation  at  the  bottom  of  the  pulses. 
Furthermore,  it  can  reduce  the  effective  width  of  the  output  pulses,  as  shown 
in  Figure  11. 


Figure  11.  Shaving  the  Wave  by  Using  a  Schottky  Diode 
3.  Results  of  Measurements 

The  output  wave  forms  of  the  control  pulse  generator  developed  at  various 
loads  are  shown  in  Figures  12  and  13.  The  wave  form  generated  on  a  pure 
resistive  load  (50Q)  is  fairly  symmetric  with  negligible  oscillation  (see 
Figure  12).  When  it  is  triggered  by  a  5  ns ,  5  V  pulse,  the  delay  time  of  the 
generator  is  5  ns.  If  a  trigger  pulse  with  a  steeper  leading  edge  is  used, 
then  the  delay  will  be  even  shorter.  When  a  40  pF  capacitor  is  connected 
parallel  to  the  output  load  (50fi) ,  the  effect  on  the  output  wave  form  is  not 
significant  (see  Figure  13). 
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Figure  12.  Wave  Form  of  the  Generator 
on  a  Pure  Resistive  Load 
(50fi);  Y:25V/div  and 
X:2  ns/div 


Figure  13.  Effect  of  Capacitive  Load 
on  the  Wave  Form  of  the 
Pulse  Generator; 

Y:25V/div  and  X:2  ns/div 


4.  Connecting  the  Pulse  Generator  to  the  ITT-F4111  Image  Intensifier 

If  the  pulse  generator  is  directly  connected  to  the  ITT-F4111  image  intensifier 
as  shown  in  Figure  7,  the  result  is  still  not  satisfactory  because  the  load  in 
this  situation  is  no  longer  a  purely  resistive  50ft  load.  Instead,  it  is  the 
impedance  between  the  photocathode  and  the  microchannel  plate  input  of  the 
image  intensifier.  As  we  mentioned  earlier,  R^  =  4  x  and  C^  =  23  pF. 

In  addition,  based  on  the  relative  optical  gain  Vj,  vs  photocathode  voltage  Uc 
curve  shown  in  Figure  2,  we  know  that  the  optical  gate  may  be  opened  to  vari¬ 
ous  extent  when  only  a  few  volts  is  applied  to  the  photocathode.  If  the  wave 
form  of  the  control  pulse  is  not  an  ideal  triangle,  then  the  oscillation  super¬ 
imposed  to  the  bottom  of  the  pulse  may  lead  to  the  continuous  opening  of  the 
optical  gate,  resulting  in  multiple  exposures.  This  must  absolutely  be 
avoided.  To  this  end,  another  pulse  shaping  and  dc  bias  circuit  is  intro¬ 
duced  between  the  pulse  generator  and  the  image  intensifier  to  result  in  the 
complete  circuit  for  a  photocathode  controlled  image  intensifier,  as  shown  in 
Figure  14. 


Figure  14.  Complete  Circuit  of  a  Photocathode  Controlled  Image  Intensifier 
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This  circuit  has  the  following  functions.  On  one  hand,  for  direct  current, 
resistors  Rx,  R2  and  R3  constitute  a  dc  voltage  divider  which  is  used  to  shift 
the  reference  voltage  (i.e.»  to  change  the  voltage  at  point  A  in  Figure  14)  of 
the  pulse  output.  Consequently,  the  reference  voltage  of  the  output  pulse  is 
+UA,  instead  of  0.  The  function  of  the  voltage  biasing  circuit  is  shown  in 
Figure  15.  Figure  15(a)  shows  that  the  oscillating  amplitude  at  the  bottom 
of  the  pulse  is  sufficient  to  cause  the  optical  gate  to  open  many  times  when 
the  reference  voltage  of  the  pulse  output  is  0V.  In  Figure  15(b),  with  the 
aid  of  the  voltage  divider,  the  reference  voltage  is  made  to  be  +UA.  As  long 
as  UA  is  sufficiently  large,  it  is  possible  to  reliably  prevent  the  opening 
of  the  optical  gate  due  to  the  oscillation  at  the  bottom  of  the  control  pulses. 
On  the  other  hand,  in  terms  of  alternating  current,  the  right  side  of  the 
coupled  capacitor  of  the  pulse  generator  (Cj^  in  Figure  14)  is  an  ac  load. 
Therefore,  the  output  pulse  is  not  affected  by  the  large  dc  resistance  between 
the  electrodes  in  the  image  intensifier.  The  nonlinear  circuit  consisting  of 
the  diode  D 3  and  R4  and  C3  is  used  to  further  shape  up  the  pulse.  D3  is  a 
Schottky  diode  with  an  extremely  low  capacitance.  The  optimal  values  of  R4 
and  C3  are  determined  experimentally. 

Another  function  of  the  voltage  biasing  circuit  is  to  adjust  the  effective 
width  of  the  output  pulse,  resulting  in  a  shorter  exposure  time.  Of  course, 
it  is  done  at  the  expense  of  some  optical  gain  because  the  effective  ampli¬ 
tude  of  the  pulse  is  also  decreased  with  decreasing  effective  pulse  width. 

The  super-high-speed  camera ,  comprised  of  a  control  pulse  generator  and  an 
image  intensifier  as  the  optical  shutter,  which  is  capable  of  capturing  faint 
instantaneous  single  events  with  a  response  in  the  nanosecond  range,  has 
already  been  used  in  the  study  of  high-speed  physical  phenomena. 


Figure  15.  Function  of  the  Voltage  Biasing  Circuit 


This  work  was  performed  at  the  Institute  of  Electrical  Measurement  of 
Darmstadt  Poly technical  University  in  West  Germany. 
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cision  Mechanics,  CAS;  manuscript  received  28  August  1984] 

[Text]  Abstract.  The  pulsed  avalanche  discharge  process 
in  an  XeCl  laser  at  pressures  up  to  10  atm  is  analyzed  and 
the  effect  of  the  delay  phase  on  discharge  formation  time 
is  discussed.  By  a  comparison  with  experimental  results 
for  a  pulsed  avalanche  discharge  in  an  XeCl  laser,  a  critical 
avalanche  track  length  of  about  1  mm  is  derived. 

Avalanche/ Self-sustained  discharges  are  extensively  used  in  excimer  lasers; 
analysis  of  the  physical  principles  of  this  type  of  discharges  at  high  gas 
pressures.  Palmer  [1]  established  a  simple  physical  model  for  estimating  the 
minimum  preionization  density  for  uniform  glow  discharges  at  high  gas  pres¬ 
sures,  in  which  the  pulse  voltage  is  applied  at  both  edges  of  the  discharge 
electrodes  as  a  stepwise  increasing  function;  shortly  thereafter  Levatter  and 
Lin  [2]  considered  the  effect  of  finite  voltage  rise  time  on  the  discharge 
process,  using  working  gas  pressures  below  3  atm.  Below  we  analyze  the  pulsed 
avalanche  discharge  process  at  pressures  up  to  10  atm  and  discuss  the  effect 
of  the  delay  phase  on  the  discharge  formation  time;  through  a  comparison  with 
experimental  results  for  a  pulse  avalanche  discharge  in  an  XeCl  excimer  laser 
we  derive  a  critical  track  length  of  about  1  mm  for  the  avalanche  discharge. 

I.  Pulse  Avalanche  Discharge  Experiments  at  High  Gas  Pressures 

In  our  experiments  we  used  a  small  discharge  apparatus  driven  by  an  X-ray 
preionization  pulse  forming  network,  whose  construction  and  voltage  measure¬ 
ment  technique  are  described  in  Ref  3. 

We  first  used  a  Rogowski  coil  to  determine  the  discharge  voltage  waveform  at 
a  total  gas  pressure  of  5  atm.  The  pulse  forming  network  consisted  of  four 
high-voltage  coaxial  cables  with  characteristic  resistances  of  50  ohms, 
connected  in  parallel  to  give  a  characteristic  resistance  of  12.5  ohms.  The 
gas  discharge  area  was  2  cm^  and  the  electrode  spacing  was  1  cm.  Discharge 
circuit  analysis  indicated  that  when  the  load  resistance  R(t)  was  less  than 
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Figure  1. 

Voltage  waveform  of  pulsed  avalanche 
discharge  in  XeCl  laser,  showing 
impedance  matching  effect  (time  scale 
50  ns  per  division;  vertical  scale 
6  kV  per  division) 


the  characteristic  resistance  Z  of  the  pulse  shaping  network,  the  voltage 
waveform  consisted  of  damped  oscillations  whose  second  peak  amplitude  U9  was 
related  to  its  first  peak  amplitude  Ui  by  the  ratio  1 


Uj  =  R—Z 
R  +  Z' 


(1) 


The  negative  sign  on  the  left  side  of  this  equation  indicates  that  U9  and  Ui 
^rev°ft  °Pposite  siSn.  Figure  1  shows  the  discharge  voltage  waveform  for  an 
Ne:Xe:HCl  ratio  of  989:10:1.  From  the  ratio  of  U2  to  Ui  we  can  determine  the 
average  value  of  R,  0.67Z  -  8.4  ohms;  this  is  rather  close  to  the  steady-state 
resistance  of  7.6  ohms  that  we  reported  in  Ref  3.  The  above  result  indicates 
that  the  line  mismatch  is  not  serious. 


Under  the  working  conditions  described  above,  we  determined  the  voltage  wave¬ 
form  for  a  total  gas  pressure  of  2-10  atm.  It  can  be  seen  in  Figure  2  that  at 
5  atm  or  less  the  voltage  waveform  agrees  with  the  analysis  in  Ref  2:  when  the 
pulse  voltage  has  risen  linearly  to  a  certain  value,  the  laser  medium  rapidly 
becomes  conductive  and  passes  over  to  a  pseudo-steady-state  discharge.  If  we 
define  the  time  elapsing  between  the  maximum  voltage  and  the  steady  state  as 
the  discharge  buildup  time  tp,  then  for  P>5  atm,  tp>20  ns,  while  when  P>5  atm 
e  voltage  waveform  changes:  there  is  a  delay  phase  between  the  time  that 
the  voltage  applied  to  the  discharge  gap  rises  rapidly  and  the  phase  during 
which  the  gas  is  conductive;  the  delay  becomes  more  pronounced  as  the  gas 
pressure  is  increased.  When  this  delay  phenomenon  occurs,  tD  also  increases 
with  increasing  gas  pressure.  Figure  2  also  shows  the  range  of  tn  for  certain 
gas  pressures;  for  clarity,  tp  is  given  only  for  total  gas  pressures  of  2,  6 
and  10  atm.  Because  the  discharge  process  is  somewhat  random,  the  values 
presented  in  Figure  2  are  a  group  of  typical  values.  It  will  be  seen  below 
that  the  experimental  data  are  obtained  by  averaging  the  values  from  5  or  more 
experimental  photographs. 


In  addition,  our  experimental  results  also  indicated  that  the  delay  phenomenon 
is  inversely  related  to  the  initial  voltage  of  the  pulsed  discharge  source,  de¬ 
creasing  when  the  source  voltage  rises.  But  if  an  excessively  high  initial 
voltage  issued,  the  laser  efficiency  may  be  extremely  low.  In  ordinary  trans¬ 
verse-discharge  excimer  lasers,  the  usual  initial  voltage  is  20-30  kV  per 
centimeter  of  gap  width  in  order  to  keep  the  laser  in  relatively  high-efficiency 
operation.  Since  we  studied  the  delay  phenomenon  under  the  conditions  described 
above,  we  kept  the  initial  discharge  voltage  constant  and  studied  the  delay 
phenomenon  under  these  constant-voltage  conditions. 
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Figure  2. 

Voltage  waveform  for  pulsed  avalanche 
discharge  in  XeCl  laser  at  10  atm 
(vertical  scale,  13  kV  per  division). 


o  100  200  300 

(ns) 


The  above  experimental  results  indicate  that  as  the  working  gas  pressure  of 
the  laser  is  increased,  the  pulsed  avalanche  discharge  process  acquires  new 
characteristics. 


II.  Analysis  of  the  Discharge  Formation  Process  at  High  Gas  Pressure 


According  to  the  general  theory  of  gas  discharges,  there  are  two  discharge 
mechanisms:  the  slow  process,  also  called  the  Towsend  mechanism,  which  does 
not  take  account  of  the  space  charge  field  produced  by  charged  particles,  and 
the  fast,  streamer  process,  also  called  the  channel  mechanism,  which  takes 
account  of  a  strong  space  charge  field  and  passes  over  directly  from  a  pri¬ 
mary  electron  avalanche  into  a  streamer  discharge;  if  this  discharge  is  not 
controlled,  it  readily  becomes  a  spark  channel.  In  order  to  suppress  the 
development  of  a  spark  channel,  the  preionization  technique  can  be  used. 
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Based  on  the  analysis  in  Ref  2,  when  the  applied  voltage  rises  linearly,  the 
electron  avalanche  head  radius  rc  at  the  initial  breakdown  is  approximately 


'c=«0Uer,  (2) 

where  £c  is  the  critical  discharge  track  length  and!e  is  the  electron  mean 

free  path.  In  order  to  obtain  a  uniform  discharge,  the  streamer  discharge 

must  develop  uniformly  in  each  part  of  the  discharge  area,  which  requires 

that  the  average  initial  electron  density  n  in  the  cathode  surface  layer 

satisfy  the  condition  n  >(r)-3.  eU 

eu 

As  the  total  gas  pressure  is  increased,  the  electron  mean  free  path  H  de¬ 
creases,  so  that  re  falls  as  the  gas  pressure  rises.  In  the  present  approx! 
mat ion,  we  ignore  the  change  in  neo  with  gas  pressure;  then  the  discharge 
buildup  time  tp  is  given  by 


where  ue  is  the  electron  drift  velocity  in  the  discharge.  According  to 
Boltzmann  coding  calculations  for  XeCl  excimer  lasers,  ue  is  a  function  of 
E/n,  such  that  when  E/n  increases  from  10“16  to  5  x  10-1°  V-cm^,  u  rises 
from  8  x  10  to  3.2  x  107  cm/sec;  the  relationship  can  be  expressed  described 
by  an  approximately  linear  formula.  In  our  experiments,  as  the  gas  pressure 
was  increased  the  number  of  particles  n  per  unit  volume  rose;  in  addition 
because  the  breakdown  voltage  increased  and  the  delay  was  lengthened ,  the ’ 
average  field  intensity  E  in  the  working  gas  before  breakdown  increased  as  a 
function  of  the  gas  pressure,  but  the  rate  of  increase  of  E  was  much  less 
than  the  rate  of  increase  of  n,  so  that  E/n  fell  as  the  gas  pressure  rose, 
further  lowering  the  electron  drift  velocity  in  the  avalanche  process. 


For  an  Ne:Xe:HCl  ratio  of  989:10:1,  we  found  the  change  in  t 
of  at  different  gas  pressures  from  Figure  2.  For  the  sake  o5 


as  a  function 
accuracy,  we 


averaged  the  results  from  five  similar  experimental  photographs.  These  re¬ 
sults  are  shown  in  Figure  3  by  the  dotted  line.  However,  when  we  calculated 
the  values  of  E/n  and  ue  for  different  gas  pressures,  we  obtained  the  inter- 
esting  result  that  the  product  of  t  and  ue  was  essentially  unchanged  at 
different  gas  pressures,  which  means  that  the  critical  track  length  was 
not  sensitive  to  changes  in  gas  pressure  (under  our  experimental  conditions 
as  escribed  above).  The  values  of  tp  calculated  for  various  gas  pressures 

from  equation  (3)  with  ?c  =  1  ™  are  shown  in  Figure  3;  the  agreement  is 
excellent . 

Figure  3. 

Discharge  formation  time  vs  total  gas 
pressure.  Solid  line  gives  theoretical 
calculated  value,  dotted  line  gives  ex¬ 
perimental  results 


V  1  3  5  7 

Total  gas  pressure 
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Ill .  Conclusions 

The  analysis  above  indicates  that  when  the  total  gas  pressure  exceeds  5  atm, 
the  delay  in  discharge  initiation  has  a  rather  large  effect  on  the  discharge 
buildup  time;  this  delay  phenomenon  is  statistical  in  nature  and  the  delay 
time  is  inversely  proportional  to  the  probability  W  that  a  primary  electron 
will  start  an  avalanche.  W  in  turn  increases  as  the  gap  overvoltage  rises. 

In  the  experiments  described  above,  in  order  to  maintain  a  constant  stored 
energy  value,  we  maintained  the  power  supply  voltage  constant.  In  order  to 
decrease  the  effect  of  this  delay  phenomenon,  the  voltage  must  be  increased 
at  high  gas  pressures  so  that  it  operates  in  a  good  overvoltage  state,  but 
this  greatly  increases  the  stored  energy  in  the  power  supply  and  decreases 
the  laser  efficiency.  In  order  to  overcome  this  conflict,  the  two-pulse 
technique  can  be  used,  in  which  a  high-overvoltage  pulse  can  be  used  to  de¬ 
crease  the  delay,  immediately  followed  by  a  lower-overvoltage  pulse  to  increase 
device  efficiency.  This  approach  has  already  been  used  successfully. 

On  the  other  hand,  the  above  comparison  of  theory  and  experimental  results 
gives  a  critical  avalanche  track  length  of  about  1  mm  for  a  pulsed  avalanche 
discharge  in  an  XeCl  excimer  laser.  This  result  is  in  agreement  with  the 
value  £c  =  1-2  mm  reported  in  Ref  2  for  mixed  gas  in  an  XeF  laser  with  a 
voltage  rise  time  of  between  10  and  20  ns  and  thus  provides  an  indirect  method 
of  determining  £c. 

Comrade  He  Qisheng  [0149  7784  3932]  took  part  in  the  experiments,  and  we 
express  our  thanks  to  him. 
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PSEUDOINVERSE  RESTORATION  OF  BLURRED  IMAGE  IN  UNIFORM  MOTION 

Chengdu  SICHUAN  UNIVERSITY  (ZIRAN  KEXUE  BAN)  [JOURNAL  OF  SICHUAN  UNIVERSITY 
(NATURAL  SCIENCE  EDITION)]  in  Chinese  No  2,  1985  pp  44-49 

[Article  by  Si  Xianyu  [5685  7359  3254],  Zhang  Guanshen  [1728  0385  3947],  Zhang 
Shaoying  [1728  1421  4481]  and  Li  Yongguo  [2621  3057  0948],  manuscript  received 
on  22  December  1983:  "A  Pseudoinverse  Restoration  Method  for  Blurred  Image  in 
Uniform  Motion"] 

[Text]  Abstract:  In  this  paper,  the  mistakes  in  the  itera¬ 
tion  formula  for  the  restoration  of  blurred  images  in  uniform 
motion  in  reference  [1]  are  identified.  The  correct  formula 
is  given.  In  addition,  a  pseudoinverse  restoration  method, 
which  adds  a  corrected  image  to  a  basic  image,  is  introduced. 

Results  obtained  from  the  computer  show  that  this  method  is 
effective. 

1.  Introduction 

In  the  image  forming  process,  due  to  various  reasons,  there  is  a  relative  motion 
between  the  recording  device  and  the  image.  For  example,  if  the  camera  moves 
relative  to  the  scene,  the  photograph  taken  will  become  blurred.  Restoration 
of  blurred  images  due  to  uniform  motion  is  one  of  the  topics  in  image  process¬ 
ing.  An  iteration  formula  for  the  restoration  of  blurred  images  in  uniform 
motion  was  discussed  in  the  literature,1  Because  of  a  mixup  of  variables, 
however ,  a  wrong  equation  was  derived.  In  this  paper,  the  iterative  restoration 
formula  was  re— derived.  In  addition,  based  on  the  concept  of  pseudoinverse 
restoration,  a  pseudoinverse  restoration  method  involving  adding  a  corrected 
image  to  the  basic  image  was  introduced.  The  feasibility  of  this  method  was 
proven  by  simulations  performed  on  a  microprocessor. 

2.  Re-derivation  of  the  Iteration  Formula 

Let  us  assume  that  the  original  image  f(x)  is  moving  uniformly  along  a  straight 
line  in  the  x-direction,  the  length  of  the  frame  is  L,  and  the  image  moves  by  a 
distance  a  within  the  exposure  time  T.  Assuming  L  =  ka  and  k  is  an  integer, 
then  the  blurred  image  becomes 

f(x-*t)di  (i) 
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A  discrete  formula  is  used  to  express  this  blurring  process.  Let  us  choose 
the  frame  length  to  be  L  image  elements  and  the  distance  of  the  image  movement 
to  be  a  image  elements ,  then 


1  - 

g(*)= — 51 /<*-«),  x  =  0,  1,  —  L-l 

a  TS 

g/(*)  =  g(*)-g(x-l),- 


(2) 

(3) 


is  used  to  represent  the  partial  difference  of  the  blurred  image.  By  substl 
tuting  equation  (2)  into  (3)  we  get 


f(x)  =  ag'(x)  +f(x-a ) 

Let  us  introduce  two  variables:  Z  and  m.  x  =  Z  +  ma,  and  m  is  the  integer 
part  of  x/a.  Z  =  0,  1,  •••  a  -  1.  The  above  equation  can  be  written  as: 

/(Z+m«)-tfg/(Z+ma)+/[Z  +  (m-l)tf]  (4) 

Let  cJ>(Z)  be  the  portion  moved  into  the  frame  from  the  outside  during  the 
exposure  period,  i.e.  (j)(Z)  =  f(Z  -  a).  Hence 

When  m  =  0,  f(Z)  =  <7g'(Z) +<£(.Z) 

n  KZ+a^-ag'CZ+a)  +/(Z)  =  og'(-Z.+  ff) +ffg'(Z) +0(Z) 

m  —  l 


By  iteration,  we  get 

W  ' 

KZ +ma)  =  a'^  g' (Z +ka') +4>(Z)  (5 

In  the  range  ka  <  x  <  (k  +  l)a,  both  sides  of  equation  (5)  are  calculated.  In 
addition,  k  equations  corresponding  to  m  =  0,  1,  •••  k  -  1  are  added  on  both 
sides  to  obtain 


*-l  k-  i  * 

5J/ck+»o-®2  2  g'(Z  +  ka~)  +  k<f>(Z) 

i.e.  "*°  *-#  . 

^(2)=4S  /(Z+m«)--|-5J  (K  -  k)g' (Z  +ka) 

Substituting  it  into  equation  (5) ,  we  get 

f(x)  =  —  VJ  KZ  +  mo)  +  ~  42  “  k)8' (Z  +  ko') 


(6) 


(7) 
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The  first  summation  on  the  right  is  also  unknown.  For  large  k’s,  it  approaches 
the  average  of  f(x).  Thus,  it  is  approximately  a  constant,  A.  Hence 

Kx)=A+a^  +  CK~k)g'a+ka)  (8) 

The  result  in  reference  [1]  is  mistaken  because  variables  m  and  k  were  mixed  up. 
This  conclusion  was  confirmed  experimentally. 

3.  The  Psuedoinverse  Restoration  Method 

The  uniform  motion  blurring  process  described  above  can  be  calculated  in  an 
equivalent  vector  space.  Let  us  assume  that  the  initial  image  vector  is  f  and 
its  sample  number  is  m,  the  blurred  image  vector  is  g  and  its  sample  number  is 
n,  and  the  displacement  is  a  =  m  -  n,  then  the  uniform  speed  blurring  process 
can  be  expressed  by  the  following  vector  matrix: 

8=Hf  (9) 


where 


f=C/(0)  /(l)  ...  /<m -l)]7 

g  =  [g(0;  g(l)  ...  g(M-l)]7 


/  1 


H  = 


.1  —  1  0  \ 

1  1  ...  •  1 

1  1  ...  1 

;.\  \  \ 

i  i  ...  i 

v — . /  / . 

a  +1 


(10) 


H  is  ann  x  m  right  shifting  cyclic  matrix.  The  restoration  problem  is  reduced 
to  solving  the  matrix  equation.  Because  m  >  n,  the  number  of  unknowns  is  larger 
than  that  equations.  Therefore ,  there  are  infinite  number  of  solutions.  Hence, 
it  is  difficult  to  restore  accurately.  Based  on  the  concept  of  a  generalized 
inverse  matrix,  the  least  modular  solution  fmi  of  the  image  vector  f  can  be 
determined. 

fru~BKHiFylg  (I1) 


Because  the  original  image  does  not  necessarily  have  the  least  amount  of  energy, 
however,  the  least  modular  solution  is  not  an  accurate  estimation  of  the 
original  image. 
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Based  on  the  concept  of  psuedoinverse  restoration, 3  a  method  to  add  a  corrected 
image  to  the  basic  image  is  introduced.  The  general  solution  to  equation  (9) 
can  be  expressed  as 


f-  Hsgq(,Im  -  HR0H)y 


(12) 


where  HR  is  a  right  inverse  matrix  of  H,  Ija  is  a  unit  vector,  and  y  is  an 
arbitrary  m  x  1  column  vector.  Let  us  denote  the  right  inverse  matrix  of  H 
whose  first  a  rows  of  elements  are  zero  to  be  HR0,  then  equation  (12)  can  be 
rewritten  as 

f-  Hsog  + ( /„  -  HR0H)y 

-/„+/»  <13) 


We  can  also  prove  that 

f  0  0  0  0 


0  0  0  0 

1 

-1  1  • 

0  -1  1 


H*°  = 


0  -1 


0 

0 

-1 

\ 

1 

0 

0 

\ 

\ 

-1 

1 

0 

\ 

\ 

\ 

0 

-1 

1 

0 

. 

\ 

\ 

\ 

• 

• 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

• 

• 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

0 

\ 

\ 

\ 

\ 

\ 

\ 

A 

0 

V _ 

•  •• 

•  •• 

0 

-1 

1 

0 

\ 

\ 

m 


(14) 


a  + 1 


a  +  1 


V. 


53 


Furthermore,  we  find 


(15) 


We  can  see  that  the  last  m  -  a  elements  are  unrelated  to  the  result  of  equa¬ 
tion  (13).  Let  y  =  ya  =  [y (1)  y(2)  •••  y(a)  00  •••  0]^.  Thus,  equation  (13) 
can  be  simplified  as 

f-  H*°y  +  ( /„'  - 

=  (16) 

The  physical  significance  of  the  above  equation  is  very  obvious.  The  first 
term  fQ  is  the  special  solution  of  the  vector  matrix  equation  (9) ,  which  is 
the  basic  image.  The  first  a  image  elements  of  the  basic  image  are  0.  This 
is  a  result  of  choosing  the  blurred  image  g  as  the  blurring  of  a  specific  part 
of  the  original  image  outside  the  frame  with  some  zero  elements.  The  second 
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term  fa  represents  the  correction  of  the  basic  image  based  on  an  estimated 
value  ya  for  the  part  outside  the  frame  when  that  part  is  not  zero.  It  is 
called  the  correction  image.  Because  the  basic  and  correction  images  are 
separated,  it  is  possible  to  estimate  the  portion  outside  the  frame  according 
to  different  criteria  to  conduct  an  interactive  restoration.  Notice  that  the 
elements  of  and  (Im  =  H^H)  only  include  0  and  +1,  it  is  possible  to  avoid 
matrix  multiplication  when  using  equation  (16)  in  the  computation.  Simple 
addition  and  subtraction  can  be  used  directly  to  drastically  reduce  the  work 
load. 

4.  Experimental  Results  and  Discussion 


This  method  was  simulated  on  a  microprocessor.  A  FORTRAN  subroutine  was  used 
to  display  the  image  as  a  gray  toned  image  given  by  a  conventional  printer. 
Different  characters  are  overlapped  to  print  out  different  tones  of  gray.  The 
results  are  shown  in  Figure  1.  Figure  (a)  shows  the  initial  image  and  the 
length  of  the  image  is  64  image  elements,  (b)  shows  the  blurred  image  created 
by  the  uniform  motion.  The  image  shift  distance  is  seven  image  elements.  It 
is  no  longer  possible  to  recognize  the  content  of  the  original  image  based  on 
this  figure,  (c)  shows  the  result  restored  by  using  the  method  introduced  in 
this  work.  . . . . . . . . 


Reproduced  from 
best  available  copy. 


Figure  1 
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As  for  estimating  the  portion  outside  the.  frame »  different  criteria  may  be 
adopted.  The  simplest  criteria  are  the  average  criterion  and  the  equal  gray¬ 
ness  criterion  in  neighboring  areas.  The  former  gives  an  estimation  based  on 
the  mean  value  of  the  blurred  image  and  the  latter  gives  an  estimation  based 
on  the  average  of  the  blurred  image  in  the  frame  and  the  adjacent  area  outside 
the  frame.  There  are  other  methods.  The  pseudoinverse  restoration  method 
based  on  adding  a  correction  image  to  the  basic  image  will  facilitate  the  use 
of  various  criteria  to  restore  the  image. 
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SCIENTISTS  AND  SCIENTIFIC  ORGANIZATIONS 


FUJIAN  MATERIALS  STRUCTURE  INSTITUTE  DESCRIBED 

Beijing  HUAXUE  TONGBAO  [CHEMISTRY]  in  Chinese  No  7 ,  18  Jul  85  pp  58-59 
[Article  by  Wu  Dingming  [0702  7844  6900]] 

[Text]  The  Fujian  Materials  Structure  Institute  of  the  Chinese  Academy  of 
Sciences  (FMSI  hereafter)  is  located  at  the  western  suburb  of  Fuzhou  on  the 
banks  of  the  Min  Jiang.  Situated  at  the  foot  of  a  hill  and  beside  a  stream, 
it  commands  magnificent  views.  It  was  established  in  1960  by  consolidating  the 
six  institutes — Technological  Physics,  Applied  Chemistry,  Electronics, 
Mathematical  Mechanics,  Automation  and  Rare  Metals — earlier  planned  and  estab¬ 
lished  by  the  Fujian  branch  of  the  Chinese  Academy  of  Sciences.  In  September 
1961,  it  was  further  consolidated  and  reorganized  into  the  Institute  of  Physi¬ 
cal  Sciences  in  accordance  with  the  then  enacted  Eight  Character  Policy.  In 
January  1962,  it  was  renamed  Huadong  Materials  Structure  Institute  and  was 
subordinate  to  the  Huadong  branch  of  the  Chinese  Academy  of  Sciences.  Its 
name  was  changed  to  the  Institute  No  7701  of  the  Office  of  Defense  Industry, 
Fujian  Province  in  1970  and  was  finally  settled  in  1973  as  the  Fujian  Materials 
Structure  Institute  of  the  Chinese  Academy  of  Sciences. 


Ourdoor  scene  of  FMSI 


Professor  Lu  Jiaxi  looking  up 
literature 


Through  20  some  years  of  development,  FMSI's  own  characteristic  research 
emphasis  has  gradually  evolved.  They  include,  within  the  area  of  structural 
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chemistry,  the  exploration  of  the  relationship  between  the  macroscopic  proper¬ 
ties  (chemical  and  physical)  and  the  crystal  and  molecular  structure  of  new 
type  compounds;  systematically  carrying  out  the  basic  and  applied  research  aim¬ 
ing  at  transition  metal  complexes  (particularly  atom  cluster  compounds) , 
biological  macromolecules  and  crystals  as  primary  targets;  properly  probing 
structure  research  methodologies  and  raising  research  standards;  development 
of  catalysts,  new  technology  crystals  and  the  applied  and  developmental 
research  on  metal  corrosion  and  its  prevention  in  light  of  the  economic  con¬ 
struction  needs;  and  emphasizing  the  combination  of  theory  and  practice, 
chemistry  and  physics  as  well  as  structure  and  property  in  order  to  build 
step-by-step  a  prominent,  integrated  research  institute  of  structural 
chemistry. 

Based  on  its  developmental  goals  and  research  responsibilities,  FMSI  has  set 
up  laboratories  in  the  areas  of  structural  chemistry,  new  technology  crystals; 
complex  catalysis,  analytical  chemistry  and  the  Shamen  Sea  Water  Corrosion 
Experimental  Station.  There  are  also  research  supporting  units  such  as 
computer  station  and  books  and  information  service.  Furthermore,  there  are 
the  small-scale  production  workshop  of  new  technology  crystals  that  have  been 
through  the  research  stage,  the  machine  shop  to  meet  the  scientific  research 
needs  and  a  general  service  system  that  serves  the  whole  institute.  The 
institute  currently  employs  more  than  570  people,  including  over  370  research 
and  technical  personnel,  in  which  13  are  senior  researchers  and  over  200  are 
research  assistants  and  engineers .  Current  director  of  the  institute  is 
Liang  Jingkui  [2733  2417  7608],  a  research  fellow. 

Lu  Jiaxi  [4151  0857  6932],  the  renowned  structural  chemist  and  the  president 
of  the  Chinese  Academy  of  Sciences,  is  the  founder  of  FMSI  and  was  its  direc¬ 
tor  until  1984  when  he  stepped  down  and  served  as  its  honorary  director. 
Professor  Lu  studied  at  the  London  University  of  England  in  his  early 
years  and  obtained  a  PhD.  From  1939  to  1945,  he  was  engaged  in  research 
at  the  California  Institute  of  Technology  and  the  University  of  California 
of  the  United  States  and  published  several  papers  on  radiation  chemistry  and 
structural  chemistry.  He  returned  in  1945  and  was  involved  in  teaching, 
scientific  research  and  the  leadership  role  in  the  teaching  body  at  the  Shamen 
and  Fuzhou  University.  He  has  trained  for  the  country  a  number  of  talents  in 
chemistry  and  physics,  many  of  whom  have  become  the  key  members  in  scientific 
research  and  teaching.  In  1960,  he  started  FMSI  and  poured  in  a  tremendous 
amount  of  effort  for  the  planning,  construction,  development  and  growth  of  the 
institute.  Now  he  is  still  in  charge  of  and  cares  about  all  aspects  of  the 
institute.  He  also  directs  the  key  research  projects  of  the  academy  such  as 
the  structural  model  study  of  the  active  center  of  nitrogen  fixation  in  nitro- 
genase  and  synthesis  of  its  analogs  and  the  synthetic  chemistry  and  structural 
chemistry  of  molybdenum,  iron  and  Mo-Fe  cluster  compounds  (particularly  the 
sulfur  containing  cluster  compounds). 

FMSI  has  all  kinds  of  modern  research  tools  and  testing  capability  of 
structural  chemistry.  There  are  large  precision  instruments  such  as  the 
PDP-11/70  computer,  FTS-20E/D-V  Fourier-transformed  Infra-red  spectrophoto¬ 
meter,  PSM-50L  ultraviolet  spectrophotometer,  PE- 306  atomic  absorption 
spectrophotometer,  G.P3.5-D  high-frequency  ionic  emission  spectrophotometer, 
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Model  1106  elemental  analyzer,  Finnigan-MAT  312  mass  spectrometer,  FT-80A 
nuclear  magnetic  resonance  spectrometer,  B-ER  420  electron  paramagnetic 
resonance  spectrometer,  CAD4-SDP34C  bicircular  diffractometer,  200mA  rotating 
anode  x-ray  machine,  D/MAX- YA  rotating  target  x-ray  machine,  AD-1  automatic 
low-light  densitometer,  DX-3A  scanning  electron  miscroscope,  CA25  4096 
channels  pulse  analyzer,  OCD  circular  dichroism  spectrometer,  and  the  home¬ 
made  Mossbauer  spectrometer,  positron  annihilation  spectrometer,  magnetic 
balance,  Model  41  differential  precision  thermobalance,  etc.  These  instru¬ 
ments  not  only  help  accomplish  the  key  projects  of  the  institute  but  also  pro¬ 
vide  services  for  the  assay  and  examination  of  various  natural  chemical  sub¬ 
stances,  crystals  and  chemical  products  as  well  as  information  for  product 
analysis,  property  clarification  and  analog  synthesis. 

Since  its  establishment,  FMSI  has  claimed  over  120  scientific  and  technologi¬ 
cal  accomplishments,  22  of  them  have  won  the  significant  scientific  and 
technological  achievement  award  of  the  National  Science  Meeting  and  of  the 
Chinese  Academy  of  Sciences.  Examples  are  the  study  on  the  simulation  of  the 
nitrogen  fixation  active  center  of  nitrogenase,  the  attempted  synthesis  of 
the  "Fuzhou  Model  II"  model  compounds,  the  improvement  of  the  catalyst  for 
the  synthesis  of  butanol  by  the  Reppe  carbonylation,  the  exploration  of  new 
materials  and  the  functional  group  model  of  nonlinear  optical  materials,  the 
study  of  yttrium  aluminate  (YAP)  laser  crystal  and  high  power  continuous 
laser,  Ti-Pb02  electrode  with  sandwiched  micro  layer  of  platinum,  introduction 
of  aldehyde  functional  groups  to  vinylon,  and  stainless  steel  equipment  and 
monitor  technology.  Many  optical  crystals  that  have  been  successfully  devel¬ 
oped  include  the  electrophoto  and  nonlinear  optical  crystal  of  ammonium 
dihydrogen  phosphate  (ADP) ,  potassium  dihydrogen  phosphate  (KDP)  and  deuterated 
potassium  dihydrogen  phosphate  (DKDP) ,  the  yttrium  aluminate  laser  crystal  with 
single  dopant  (neodymium)  or  double  dopants  (neodymium ,  chromium)  (Nd,  Nd+Cr, 
YAP) ,  the  x-ray  diffraction  crystals  such  as  potassium  hydrogen  phthalate 
(KHP) ,  Pentaerythritol  (PET) ,  lead  stearate  (LS)  and  lead  cerotate  (LCer) ,  and 
the  crystal  growth  of  the  protein  from  Trichosanthes  kirilowii  and  its 
crystal  cell  coordinates  and  low  resolution  structure  determination  (in 
collaboration  with  the  Shanghai  Organic  Chemistry  Institute  and  the  Institute 
of  Biophysics).  There  are  also  the  new  nonlinear  optical  crystal,  the  low- 
temperature  phase-polarized  barium  borate,  that  has  been  discovered  and  suc¬ 
cessfully  grown  after  many  years  of  efforts.  In  recent  evaluation  by  the 
academy,  this  is  regarded  as  a  successful  development  in  the  research  of 
new  technology  crystals  solely  by  our  own  scientific  and  technical  personnel 
through  combining  theory  with  practice  and  physics  with  chemistry  and  one  that 
has  left  our  mark. for  the  first  time  on  the  frontier  of  material  science.  It 
is  believed  to  have  significant  impact  on  the  advancement  of  laser  technology. 

Besides  bringing  the  key  research  projects  of  the  state  and  the  academy  to 
completion,  FMSI  also  actively  exploits  its  technical  expertise  to  serve  the 
application  needs  of  Fujian  provincial  government  and  other  production 
departments.  Several  technical  accomplishments  have  been  achieved  and  related 
production  departments  have  derived  varied  degrees  of  economic  benefits  from 
them.  Commendations  and  awards  from  the  provincial  government  and  Fuzhou  city 
as  well  as  many  other  related  departments  have  been  given  us  for  this  feat. 
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FMSI  has  been  authorized  by  the  Academic  Degrees  Committee  of  the  State  Council 
to  accept  and  train  doctoral  and  master  degree  students  and  to  award  degrees  in 
physical  chemistry  and  solid-state  physics .  At  the  moment ,  there  are  20  mas¬ 
ters  degree  candidates  and  2  doctoral  candidates. 

FMSI  has  also  actively  pursued  academic  exchange  with  foreign  countries  in  many 
forms  such  as  tour  of  investigation  in  foreign  country,  engaging  in  advanced 
studies,  attending  international  academic  meetings,  studying  for  degrees, 
collaborative  research,  short  courses  of  technical  training,  inviting  foreign 
scholars  and  experts  to  visit,  give  lectures,  work  for  short  term  or  set  up 
and  test  instruments  and  equipment.  In  recent  years,  we  have  sent  over  50 
times  all  sorts  of  scientific  and  technical  personnel  abroad  and  have  received 
on  more  than  130  occasions  visiting  experts,  which  have  positively  helped  in 
maintaining  a  lively  scientific  atmosphere  and  raising  the  academic  standard 
within  the  institute.  For  example,  Professor  Lin  Weizhen  [2651  1983  2823], 
a  physical  chemist  from  the  University  of  British  Columbia,  Canada,  was  first 
invited  in  1976  to  teach  and  work  for  a  short  term.  In  1979,  it  was  approved 
by  the  Academy  of  Sciences  to  offer  him  the  honorary  research  professorship  at 
the  institute.  Mr  Lin  has  made  significant  contributions  on  the  aspects  of 
research  guidance  and  classroom  lectures  during  his  many  short  visits  to  FMSI 
in  the  past  years.  With  the  sponsorship  of  the  academy,  he  gave  lecture 
series  at  FMSI  on  "Basic  Quantum  Chemistry  and  the  Ligand  Field  Theory" and 
"The  Irreducible  Tensor  Method"  in  1978  and  1979  that  greatly  benefited  our 
institute  and  our  sister  organizations  in  the  training  of  research  and 
teaching  staff. 

In  order  to  present  new  advancements,  new  achievements  and  new  directions  in 
each  research  area  and  to  exchange  scientific  results  and  experience,  FMSI 
started  a  new  journal  STRUCTURAL  CHEMISTRY  in  1982  on  the  foundation  of  its 
sporadically  published  FMSI  NEWSLETTER.  It  includes  as  special  issue  the 
English  edition  Crystal  Structure  Report  and  is  distributed  publicly.  Also, 
many  monographs  and  collected  papers  on  structural  chemistry  have  been  edited, 
translated  and  published  in  recent  years. 

Currently,  FMSI  is  actively  implementing  the  guiding  ideology  that  scientific 
research  be  geared  toward  the  needs  of  economic  construction  and  carrying  out 
the  policy  of  "greatly  strengthening  applied  research,  actively  and  selective¬ 
ly  participating  in  developmental  work  and  paying  continued  attention  to 
basic  research."  We  have  enthusiastically  carried  out  reform  in  order  to  make 
new  contributions  in.  answering  the  needs  of  the  four  modernizations  drive  and 
producing  more  achievements  and  talents. 

12922/9365 
CSO:  4008/1086 
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SCIENTISTS  AND  SCIENTIFIC  ORGANIZATIONS 


FIFTY  YEARS  OF  ISOTOPE  CHEMISTRY  RESEARCH  RECOUNTED 

Beijing  HUAXUE  TONGBAO  [CHEMISTRY]  in  Chinese  No  7 ,  18  Jul  85  pp  51-57 

[Article  by  Zhang  Qinglian  [1728  7230  5571]  (Footnote)  (About  the  author: 
Professor  Zhang  Qinglian  was  born  at  Changshu,  Jiangsu  in  1908.  He  studied 
at  the  graduate  school  of  Qinghua  University  from  1931  to  1934  and  obtained 
his  PhD  from  Berlin  University  of  Germany  in  1936.  The  positions  he  has  held 
include  chemistry  professor  at  the  Southwestern  Consolidated  University  and 
Qinghua  University;  professor  and  chairman  of  the  Chemistry  Department  of 
Beijing  University;  deputy  director  of  the  Chemistry  Department,  Chinese 
Academy  of  Sciences;  head  of  Stable  Isotope  Division,  State  Scientific  and 
Technological  Commission;  member  of  the  standing  committee,  the  Chinese 
Chemical  Society;  chairman  of  the  Chinese  Mass  Spectroscopy  Society  and 
member  of  the  Committee  on  Atomic  Weight  and  Isotope  Abundance  of  IUPAC. 
Professor  Zhang  specializes  in  inorganic  chemistry  and  has  long  been  involved 
in  the  teaching  and  research  in  this  area.  Since  1953,  he  has  been  working  on 
isotope  chemistry  and  heavy  water  research  and  taken  part  in  the  research  and 
development  of  various  light  isotopes)] 

[Text]  The  need  of  new  nuclear  materials  for  the  modern  science  and  technology 
of  nuclear  energy  has  mandated  the  probe  of  isotopes  of  each  element  as  a 
further  step  beyond  the  study  of  chemical  elements.  As  a  result,  the  border¬ 
line  discipline  of  isotope  chemistry  is  advancing  rapidly.  The  isotope 
chemistry  studies  the  chemical  effect  that  is  linked  to  the  change  of  electron 
energy  levels  due  to  the  difference  in  number  of  neutron  within  the  nucleus  of 
the  same  element.  The  area  of  utmost  importance  in  this  field  is  the  separa¬ 
tion  of  isotopes.  It  has  pushed  forward  the  study  of  such  chemical  effect  as 
exchange,  equilibrium  and  kinetics  as  well  as  the  applications  of  analytical 
and  tracer  methodologies. 

Since  1935 ,  the  author  has  corroborated  with  nearly  50  colleagues  on  several 
research  projects.  Brief  summaries  of  works  at  each  stage  can  be  found  in 
many  reviews.11  >17»21  »22  >38  >41  »42  »61 

On  this  occasion  of  the  50th  anniversary,  though  the  work  continues,  it  seems 
appropriate  to  sum  up  my  experience.  What  follows  is  my  attempt  to  outline 
the  four  areas  of  my  work. 


61 


I.  The  Physical  Properties  of  Isotopic  Compounds 

1.1  The  Density  and  Heat  Expansion  of  Heavy  Water  at  25°C 

The  density  of  heavy  water  at  25°C  is  an  important  physical  constant  in 
isotope  chemistry,  the  accurate  measurement  of  which  has  been  continued  for 
half  a  century.  Only  after  the  hydrogen  and  oxygen  isotopes  in  water  are 
accurately  measured  by  mass  spectrometry  can  the  density  of  heavy  water  be 
accurately  determined.  In  the  1970's,  determinations  to  seven  significant 
figures  started  to  appear.  Based  on  the  density  maximum  of  standard  mean 
ocean  water  (SMOW),  namely  999.975  kg/m3,  the  author  and  co-workers  had  made 
two  accurate  determinations  on  three  heavy  water  samples.  56  *6  7  The  four  most 
recent  determinations  of  the  density  at  25 °C  of  the  heavy  water  that  have 
the  same  composition  of  oxygen  isotopes  as  in  SMOW  are  listed  in  Table  1. 

The  results  from  three  labs  agree  well  with  one  another  and  show  that  the 
most  reliable  data  so  far  for  the  density  of  heavy  water  at  25°C  is  the 
following: 

dD2°SMow  =  1104.467(5)  kg/m3 

In  1941 ,  we  measured  the  heat  expansion  of  heavy  water  and  extended  the 
temperature  range  to  50°C,1S  which  was  later  further  extended  to  100°C  in 
1949. 26  ,27 , 31  Ajuong  the  internationally  accepted  data  of  heat  expansion, 
so  far  the  best  ones  in  the  range  of  85°  to  100°C  are  those  determined  by 
the  author  and  co-workers  (G.S.  Kell,  1977). 


Table  1. 


Year 

Authors 

Method 

Sample 

number 

kg/m3 

1975 

Babeliowsky,  Brulmans 

Large  magnetic  float 

2 

1104.475(3) 

1975 

Ceccaldi,  Girard, 
Menache,  Riedinger 

Archimedes  principles 

2 

465(3) 

1983 

Nan  Chunbo ,  Qian  Qiuyu , 
Zhang  Qinglian 

Precise  pycnometry 

2 

467(9) 

1984 

Qian  Qiuyu,  Zhang 
Qinglian 

II  If 

1 

466(9) 

1.2  The  Temperature  and  Freezing  Point  Depression  at  the  Density  Maximum  of 
Heavy  Water 

By  using  a  16ml  quartz  pycnometer,  the  temperature  at  the  density  maximum  of 
heavy  water  was  determined  in  1941  as  11.21±0.05°C. 14  Again  in  1948,  we 
determined,  using  heavy  water  solutions  of  different  acetone  concentrations, 
the  molar  freezing  point  depression  constant  of  heavy  water  as  2.00±0.01°C. 19 
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1.3  The  Vapor. Pressure,  Heat  of  Vaporization,  Critical  Constant  and  State 
Diagram  of  Heavy  Water 

Using  a  quartz  differential  barometer  and  a  paraffin  oil  bath,  we  determined 
in  1936  the  vapor  pressure  difference  of  H2O  and  D2O  in  the  temperature  range 
of  21.8°  to  232. 6 °C  and  observed  that  the  largest  difference  of  83  torr 
occurred  at  170 °C  and  that  their  vapor  pressure  was  the  same  at  225°C  and  the 
D20’s  became  greater  than  the  ^O's  beyond  this  temperature. 3  This  verifies 
the  phenomenon  that  the  critical  temperature  of  D2O  is  lower  than  that  of  H2O 
as  determined  by  the  quartz  capillary  tube  method  in  1935.  The  critical 
temperature  tc  of  D20  and  H20  mixtures  follows  the  empirical  formula  below  •  ; 

t  -■  374.2  -  2.7n 
c 

where  n  is  the  mole  fraction  of  D20. 

Comparative  determinations  of  the  critical  pressure  pc  and  the  critical 
density  dc  have  also  been  made  and  the  data  are  listed  in  Table  2. 


Table  2. 


Compound 

tc,  °C 

pc ,  atm 

dc,  g/cm3 

H20 

374.2 

218.5 

0.325 

D20 

371.5 

218.6 

0.363 

Within  the  temperature  range  of  approximately  360°  to  374. 2°C,  40  density 
readings  were  obtained  that  enabled  us  to  plot  the  density- temperature  diagram 
of  H20  and  D2O. 

Take  the  vapor  pressure  data  of  H?0  and  D20  within  the  temperature  range  of 
3.8°  to  230°C  and  plot  lg  (PH2O/PD2O)  versus  1/T,  the  difference  in  heat  of 
vaporization  can  be  calculated  from  the  slope  of  the  curve.  Some  of  the  data 
are  listed  in  Tdble  3. 


Table  3. 


H20  Heat  of 

D2O  Heat  of 

Ratio  of  heat  of 

vaporization 

vaporization 

vaporization 

°C 

cal/mole 

cal/mole 

D2O/H2O 

3.8 

10703 

11109 

1.038 

50 

10255 

10517 

1.026 

100 

9719 

9919 

1.021 

150 

9101 

9257 

1.017 

200 

8353 

8465 

1.013 

230 

7808 

7894 

1.011 
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The  difference  in  heat  of  vaporization  (D20-H20)  reduces  from  406  cal/mole  at 
3.8°C  to  86  cal/mole  at  230°,  thus  disproves  what  was  believed  in  the  litera¬ 
ture  at  that  time  that  the  difference  was  constant  between  20°  and  90°C. 

1.4  Solubility  in  Heavy  Water  and  the  Distribution  Constant  of  Solute  Between 
Two  Phases 

In  1948/49,  we  determined  and  compared  the  solubility  of  various  potassium 
salts  in  light  and  heavy  water  at  25°C.  The  results  are  listed  in  Table  4, 
in  which  S  =  moles  of  salt/55. 51  moles  of  water. 23 »24 »25  It  can  be  seen  that 
the  smallest  difference  in  solubility  is  only  5  percent  whereas  the  largest  is 
25  percent.  Also,  the  percentage  difference  of  potassium  halide  KX  solubility 
increases  with  increasing  radius  of  halide  ion.  KX03  also  behave  similarly. 


Table  4. 


Potassium 

salt 

sd2o/,sh2o 

Potassium 

salt 

sd2o/:sh2o 

KCiOtf 

0.947 

KBr03 

0.895 

KC103  . 

0.935 

KBr 

0.888 

K2Ci04 

0.932 

KI 

0.884 

KC1 

0.906  ... 

kio3 

0.833 

KRe04 

0.900 

k2so4 

0.795 

K3Fe(CN) 0 

0.897 

K2Cr207 

0.745 

Because  the  solubility  ratios  of  K2Cr207  and  Cdl2  deviate  substantially  from 
unity,  their  solubilities  Sn32»34  in  a  mixture  of  light  and  heavy  water  at 
25 °C  have  been  determined  as  follows: 

K2Cr207  Sn  =  0.5024  -  0.1467n  +  0.017n2 
Cdl2  Sn  =  2.347  -  0.7046n  +  0.068n2 

where  n  is  the  mole  fraction  of  D20.  The  phase  diagram  of  the  K2Cr207- 
H20-D20  three-component  bystem  at  25°C  has  also  been  drawn.  Similar  study 
has  also  been  done  on  NaCl12: 

NaCl  Sn  =  6.145  -  0.334n 

In  1950,  the  solubility  of  thallous  nitrate  in  heavy  water  in  the  range  of  5° 
to  50°C  were  measured30  and  the  results  disprove  the  occurrence  of  phase 
change  between  25°  and  28.4°C  as  recorded  in  the  literature. 

When  a  substance  is  dissolved  in  two  immiscible  phases ,  there  exists  a  dis¬ 
tribution  constant.  A  typical  example  is  that  the  distribution  constant  of 
iodine  in  H20  and  CC14  at  25°C  is  1:85.4,  the  concentrations  in  both  phases 
being  presented  in  grams  of  iodine/liter  of  solution.  In  1944,  the  author 
used  semi-micro  method  to  determine  the  distribution  constant  of  iodine  in 
D20  and  CC14  as* 1:103. 18 
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1.5  Thermodynamic  Properties  of  Heavy  Water  Hydrates 

By  treating  anhydrous  manganese  sulfate i  zinc  sulfate  and  cobalt  sulfate 
with  slight  excess  of  heavy  water  with  99.9  atom  percent  D,  three  newghe|vy72 
water  hydrates  were  prepared:  MnS04*4D20,  ZnS.04*7D20  and  CoS04*7D20.  * 

The  1H-NMR  was  used  to  determine  the  isotope  grade  of  the  products.  The 
differential,  the rmal/thermo gravimetric  analysis  showed  that  the  dehydration 
initiation  temperature  at  each  stage  generally  was  about  6-7° C  higher  than 
the  corresponding1 light  water  counterpart  and  in  certain  stages  the  activa¬ 
tion  energy  of  dehydration  was  higher  by  about  2-5  kcal/mole.  All  these 
support  that  the  heavy  water  hydrates  are  more  stable  than  their  light  water 
counterparts . 

1.6  The  Vapor  Pressure,  Heat  of  Vaporization  and  Boiling  Point  of  Partial 
Heavy  Water 


In  1936,  we  evaporated  at  100°C  the  water  containing  0.0347  percent  HDO  and 
measured  the  total  enrichment  of  HDO  and  H2180  in  the  residual  water  by  the 
float  method.  Then  the  H2180  content  was  determined  by  H2S  balance  method 
and  the  difference  was  HDO  content.  By  the  Rayleigh  evaporation  equation, 
the  ratio  of  vapor  pressure  HD0/H2Q  at  100°C  was  calculated:  HDO/  H2O  —  0.9734. 
From  this,  it  can  be  derived  that  the  boiling  point  of  HDO  is  100. 7°C.  Using 
the  known  data  for  lower  temperatures ,  we  derived  the  empirical  formula  for 
11°  to  100°C  range  as  follows: 


phdo/ph2o  1,16e 


/CRT) 


Hence,  the  average  difference  of  heat  of  vaporization  between  HDO  and  H20  was 
calculated  as  130  cal/mole. 4 


1.7  The  Vapor  Pressure,  Heat  of  Vaporization,  Boiling  Point,  Heat  Expansion 
and  IR  Spectrum  of  Heavy  Oxygen  Water 

In  the  above  mentioned  experiments,  we  also  measured  the  vapor  pressure  of 
H2180  at  65°,  80°  and  100°C.  Combining  this  with  the  known  data  for  lower 
temperatures,  the  following  can  be  obtained: 

„  .  -13/<RT) 

PH2180/pH20  ~  1,013 

This  is  also  applicable  from  11°  to  100°C.  It  can  bd  seen  that  the  average 
heat  of  vaporization  ofH2180  is  lower  than  that  of  H20  by  13  cal/mole. 

Table  5  lists  some  physical  constants  of  the  isotopic  species  of  water. 

The  vapor  pressure  data  of  D20,  HDO  and  H2180  measured  by  the  author  and  co¬ 
workers  have  been  applied  in  the  technology  of  heavy  water  and  heavy  oxygen 
water  productions.  '  ;  ’ ; 


Table  5. 


Isotope  substituted  water 

H20 

HD0 

D20 

H21h0 

Freezing  point,  °C 

0 

2.11 

3.80 

0.3 

Boiling  point,  °C 

100 

100. 7o 

101.40 

1001.1 3 

Vapor  pressure  at  20 °C,  Torr 

17.535 

16. 27 

15. 24 

17. 37 

Vapor  pressure  at  100°C,  Torr 

760 

739.7 

722.2 

756.5 

Heat  of  evaporation  at  100°C, 
cal  mol 

9719 

9849 

9919 

9732 

Besides,  in  1964,  we  took  a  water  sample  with  32.47  atom  percent  180  and 
treated  it  with  hydrogen  containing  0.003  atom  percent  D  in  the  presence  of 
Pt/C  catalyst  to  normalize  the  hydrogen  in  the  sample  to  150ppmD.  The  heat 
expansions  of  this  H7  80  water  sample  were  measured  by  pycnometry,  which  were 
extrapolated  to  100  percent  H2180.  Using  the  known  specific  gravity  value  of 
heavy  oxygsu  water  at  25  C,  an  empirical  formula  for  the  heat  expansion  of 
H2  80  valid  at  the  range  of  5°  to  90°C  was  derived45: 

d^  =  1.11254  +  4.16  x  10" 6 t  -  9.6  x  10-9t2 


Using  a  Nicolet  7199B  FTIR  spectrophotometer,  the  IR  absorption  peaks  of 
H2  0  and  D2180  liquids  were  observed  for  the  first  time  and  were  listed 
below  in  cm-1: 


h18o 

D180 

h218o  d218o 

3394 

2488 

1639  1202 

Further  measuring  the  IR  absorption  peaks  of  a  1:1  mixture  gave  the  V2  of 
HD180  as  1452cm-1. 


1.8  The  Boiling  Point  and  Density  of  Heavy  Ethanol 

In  1947,  we  prepared  C2H5OD  by  hydrolysis  of  sodium  ethoxide  and  aluminum 
ethoxide  with  the  heavy  water  containing  99.5  atom  percent  D  and  measured  its 
properties.  Their  comparisons  with  those  of  C2H50H  are  listed  in  Table  6. 20 


Table  6. 


Isotopic 

Density  at 

compound 

25 °C,  g/ml 

C2H5OH 

78.;32 

0.78506 

C2H5OD 

78.8 

0.801 
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I. 9  The  Shifts  in  IR  Spectra  of  Lithium  Isotopic  Compounds 

The  perchlorates  of  lithium  isotopes,  6Lio . 90 57Lio . 095CIO4  and  7Lii . 0 q 0CIO4 
and  their  crown-4,  crown-5,  crown-6,  arid  crown-8  complexes  were  prepared. 

Their  IR  absorption  spectra  from  lOO-SOOcm”1  were  studied  with  a  Nicolet 
7199B  FTIR  spectrophotometer.  Among. the  16  pairs  of  absorption  peaks  wit 
isotope  shift,  the  differences  in  wave  number  of  10  pairs  were  4.0-6. 8  per¬ 
cent,  which  is  roughly  equal  to  the  relative  difference  of  the  square  roots 
of  the  mast  of  6L1~0  and  7U-0,  i.e. .  5.5  percent.62  From  the  mid-range 
IR  spectra  of  the  crown-5  and  crown-6  complexes  of  Li  and  Lx  salts  that 
contain  water  molecules,  the  water  absorption  peaks  of  lithium  isotope 
shift  by  as  much  as  9cm"1  in  the  region  of  3436-3599cm“  . 

II.  The  Kinetic  Effect  and  Isotope  Separation  of  Isotopic  Compounds 
2.1  The  Reaction  Rates  in  Heavy  Water 

Three  reactions  were  studied  in  heavy  water.  We  designed  a  semimicro  method 
in  which  1  ml  of  reaction  solution  was  put  in  a  tiny  crucible  inside  a  con¬ 
stant  temperature  water  bath.  Heavy  water  of  different  concentrations  were 
used  and  the  data  were  extrapolated  to  100  percent  D20  medium.  At  20  C,  the 
reaction  rate  of  D202  and  V  in  D20  was  observed  to  be  only  60  percent  of 
that  of  the  corresponding  H202  and  I"  in  H20. 

Next,  we  studied  the  reaction  rate  of  D202  and  Mn04  in  acidic  D20  at  20  C. 

In  D20 ,  the  autocatalysis  of  *fri2+  ion  in  the  presence  of  D30  is  faster  than 
the  corresponding  reaction  in  H20  in  the  presence  of  H30+  while  the  rate  that 
D202  and  Mn04  react  is  only  15  percent  that  of  the  H202  and  Mn04  reaction  in 

H20. 16 

Also  the  Landolt  reaction  was  studied  at  25°C.  The  rates  that  SO |  ^reduces 
10 3  in  H20  and  D20  were  compared.  Interestingly,  it  was  observed  that  the 
rate  constant  ratio  k^/kg  was  1.2  at  0.01  molar  of  hydrogen  ion  and  was 
increased  to  3.7  at  0.5  molar. 

2.1  The  Electrolytic  Coefficient  of  Protium  and  Deuterium 

The  classical  method  of  producing  heavy  water  is  by  electrolysis,  which  was 
later  used  only  for  high  concentration  solutions.  We  studied  the  change  of 
electrolytic  coefficient  a  at  high  current  density  in  1949,  a  being  the  ■ 
relative  evolution  rate  of  protium.  A  1  percent  solution  of  D20  was  examined 
at  room  temperature.  It  was  made  basic  with  NaOH.  Using  a  platinum  wire 
with  a  diameter  of  1  mm  as  cathode,  the  experiment  was  carried  out^from  1  to 
13  amp /cm2.  It  was  found  that  a  =  6.5  ±  O^under  10  amp /cm  and  dropped 
noticeably  in  the  range  of  10  to  13  amp/cm  . 

In  1982 ,  the  fabrication  of  the  electrolytic  cell  was  improved  and  precise 
chromatographic  method  was  used  to  measure  deuterium  so  that  the  precision  ot 
a  measurement  reached  ±0.1.  An  electrolysis  experiment  was  done  with  a 
1  percent  D20  solution,  made  basic  by  K2C03.  A  regular  steel  plate  of  0.22 
percent  carbon  content  was  used  as  cathode  and  the  current  density  from  22  to 
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421  mA/cm2  was  applied.  The  a  value  was  observed  to  increase  from  7.8  to 
10.1  (at  approximately  250  mA/cm2)  and  then  dropped  to  9.5. 53 

Later,  the  a  values  of  the  high  concentration  heavy  water  were  measured. 

Under  similar  condition  but  with  197  mA/cm2  and  constant  temperature  of 
23.5°C,  the  a  was  constant  at  9.1  ±  0.1  when  the  electrolytic  solution 
contained  91.8  to  99.5  atom  percent  D.64  This  result  may  have  implications 
in  heavy  water  production  technology. 

2.3  The  Overall  Distribution  Constant  of  Deuterium  in  the  Girdler-Sulfide 
Process 

The  overall  distribution  constant  of  deuterium  3  in  the  GS  process,  a  popular 
method  for  heavy  water  production,  was  determined  in  1965.  A  teflon-sealed 
stainless  steel  exchange  chamber  of  500  ml  capacity  and  sampling  chambers  for 
both  gas  and  liquid  phases  were  fabricated.  The  heavy  water  containing  1.5 
atom  percent  D  and  the  hydrogen  sulfide  gas  of  98.6  percent  purity  were 
reacted  at  20  atm  for  5  hours  to  establish  equilibrium  under  various  tempera¬ 
tures.  The  H2S  saturated  aqueous  phase  and  the  water  vapor  saturated  H2S  gas 
phase  were  sampled  and  converted  into  pure  water,  which  were  then  analyzed 
with  precise  density  measurements.  The  experimentally-determined  6  values 
were  accurate  to  ±0.02  and  were  in  good  agreement  with  the  generally 
accepted  calculated  values  as  shown  in  Table  7. 49 


Table  7. 


Temperature 
pressure , 

30 

30 

35 

35 

130 

135 

kg /cm  2 

18 

20 

18 

20 

20 

20 

6,  obs. 

2.292 

2.33s 

2.27z 

2.260 

1.65o 

1.586 

B»  calc. 

2. 293 

2 . 28g 

2.25s 

2.254 

1.623 

1 .584 

2.4  Attempts  To  Enrich  Heavy  Water  by  Fractional  Crystallization 

In  1936,  a  copper  tube  of  20  mm  diameter  with  sealed  bottom  end  was  fashioned 
in  such  a  way  that  the  bottom  section,  15  cm  in  length,  was  bent  at  135°  angle 
and  was  loaded  with  dry  ice.  The  tube  was  spun  in  0°C  water  so  that  ice  formed 
was  quickly  pulled  down.  But  no  enrichment  of  heavy  water  in  ice  was  observed 
even  at  a  sedimentation  speed  of  1.5  m/sec.  In  fact,  we  were  limited  by 
the  precision  of  the  flotation  method  of  analysis.  Recently,  deuterium 
enrichment  in  ice  was  observed  with  a  =  1.02 

III.  Analysis  Methods  for  Stable  Isotopes  arid  the  Preparation  of  Isotope 
Reference  Samples 

3.1  The  Float  Method 

The  heavy  water  concentration  of  a  water  sample,  expressed  in  terms  of  atom 
percent  D,  can  be  measured  by  suspending  a  hollow  borosilicate  glass 
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or  quartz  float,  about  15mm  high,  in  the  heavy  water  of  certain  concentration 
at  a  specific  temperature.  It  is  a  better  way  to  measure  the  rising  and  fall¬ 
ing  velocity  of  the  float  beyond  and  below  the  flotation  temperature  to  and 
then  intrapolates  to  to.  By  using  the  density  data  of  regular  water  and 
heavy  water  at  various  temperatures  and  the  expansion  coefficient  of  the 
float,  the  heavy  water  concentration  in  a  water  sample  can  be  calculated. 

In  1956,  the  velocity-temperature  relationship  of  the  float  was  studied  and 
it  was  found  that  linear  relationship  only  existed  within  certain  temperature 
range.33  Four  types  of  float  were  later  studied  and  compared  and  the  spindle- 
shaped  float  was  found  to  be  most  satisfactory. 40  With  temperature  kept  con¬ 
stant  to  within  ±0.0005°C,  the  accuracy  of  the  measurements  reached  ±0.2  Y, 
where  y  =  mg/1.  At  25°C,  y/0.1074  =  ppmD  and  y/0.1123  =  ppm  0. 

When  analyzing  high  concentration  heavy  water,  samples  are  degassed  and 
re— dissolving  of  air  should  be  avoided  during  measurement.  Therefore,  a 
flotation  tube  with  a  ground  glass  stopper  and  a  capillary  tubing  within  has 
been  designed.57  This  is  suitable  for  the  testing  of  heavy  water  products. 
After  K2C03  is  added,  a  product  with  about  99.8  atom  percent  D  was  subjected 
to  electrolysis,  using  iron  cathode  (a  =  9.1)  and  gold  anode  (a  =  1.005  for 
160/180).  When  1/4  of  the  liquid  remained,  D  was  enriched  to  99.999  atom 
percent.  This  enriched  sample  was  distilled  out  and  compared  with  the  prod¬ 
uct  by  the  new  type  flotation  tube.  The  accuracy  reached  ±0.01  atom  percent 
D.  This  result  was  checked  in  four  labs  and  also  checked  by  mass  spectro- 
metric  method.  They  all  agreed  within  the  error  range.  Taking  into  account 
the  errors  in  heavy  water  density  data,  the  accuracy  was  determined  to  be 
±0.02  atom  percent  D.  In  1964,  when  mass  spectrometers  were  not  readily 
available,  this  was  adopted  by  the  Ministry  of  Chemical  Industry  as  a  temporary 
standard  method. 

In  order  to  extend  the  analysis  to  a  wider  range  of  concentrations  and  to 
shorten  measuring  time  to  an  hour,  a  temperature-pressure  flotation  method 
was  devised58  and  items  measured  are  listed  in  Table  8.  When  checked  by  the 
specific  gravity  bottle  method,  it  was  found  that  the  mean  deviation  was 
±0.0005  atom  percent  D  for  heavy  water  samples  of  0.11  to  0.83  atom  percent  D. 


Table  8. 


Item 

Reference  water 

Sample  water 

Thermostat  temperature,  Beckmann,  °C 

ti 

Pressure  coefficient  at  flotation 
temperature,  °C/torr 

Ki 

k2 

Flotation  pressure ,  deviation 
from  760  torr 

Pi 

P  2 

Normal  flotation  temperature, . °C 

3.2  The  Crystal  Flotation  Method 

We  studied  the  LiF  crystal  flotation  method  for  measuring  the  abundance  of 
lithium  isotopes.  The  redistilled  bromoform  was  mixed  with  n-pentanol  and 
n-hexano 1  to  get  three  mixtures  of  different  compositions  that  were  suitable 
for  measuring  different  abundances.  Their  densities  varied  from  2.63  to  2  65 
g/cm  at  temperatures  from  25.7°  to  28.8°C  and  their  stabilities  were 
examined. 

A  furnace  was  made  with  platinum  wire  that  could  be  heated  to  980°C.  A  plati¬ 
num  crucible  with  a  capacity  of  either  1  or  5  ml  was  used  to  hold  the  molten 
LiF,  from  the  surface  of  which  single  crystals  were  produced  by  using  a 
platinum  needle  1  mm  in  diameter,  which  was  soldered  to  a  platinum  tubing 
3  mm  in  diameter.  Outside  the  tubing  was  another  soldered  platinum  tubing 
6  mm  in  diameter.  Inside  the  platinum  tubing  was  a  small  glass  tubing 
through  which  the  cooling  water  flowed  to  the  bottom  and  was  directed  to  exit 
at  the  top  to  provide  cooling  for  the  needle.  The  tubing  was  fitted  with  a 
sliding  mechanism.  ,  . 

The  single  crystals  were  dome-shaped  with  2  mm  in  diameter.  The  flotation  of 
the  crystal  was  actually  measured  by  the  temperature  convergence  method.  By 
using  the  Beckmann  thermometer,  temperatures  were  measured  with  an  accuracy  of 
±0.002  C  and  ±0.004°C  for  repeated  measurements.  A  depleted  waste  product 
from  Russia  was  determined  to  contain  3.13±0.09  atom  percent  6Li.  The  result 
was  confirmed  by  the  mass  spectrometric  determination,  which  gave  3.10±0.05 
atom  percent  faLi. 51 

3.3  The  Gradient  Tube  Method 

In  1960,  we  improved  the  gradient  tube  method  by  using  new  media, 
p-dimethylbenzene  and  bromobenzene.  p-Dimethylbenzene  was  carefully  poured  on 
top  of  the  bromobenzene  layer  inside  the  gradient  tube  and  was  allowed  to 
diffuse  for  50-60  hours  at  25°C  to  form  density  gradient.  Tiny  1  yl  drops  of 
heavy  water,  whose  concentrations  were  known,  were  suspended  in  the  gradient. 

By  plotting  the  relative  height  vs  D  concentration,  standard  curves  were 
obtained  for  two  concentration  ranges  of  0-1.2  and  0-12.  The  D  concentrations 
of  testing  samples  were  then  determined  from  the  curve  by  their  heights.40 
The  accuracy  is  about  the  same  as  that  of  the  falling  drop  method. 

3.4  The  Falling  Drop  Method 

In  1957,  we  measured  the  density  of  pure  o-f luorotoluene  from  14°  to  30°C  and 
designed  a  horizontal  liquid  transfer  tube  to  produce  drops  whose  volume 
variations  were  within  0.3  percent.  The  temperature  of  the  thermostat  was 
controlled  to  within  ±0.0005°C.  The  accuracy  of  heavy  water  analysis  by  this 
method  was  improved..  When  a  one  atom  percent  D  heavy  water  sample  was 
analyzed,  the  relative  accuracy  was  improved  to  ±0.3  percent  from  ±1  percent 
attained  at  that  time.35*36 

Two  years  later,  taking  advantage  of  the  property  that  organic  liquids  have 
far  greater  heat  expansion  than  water,  we  used  a  mixture  of 
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a-chloronaphthalene  and  kerosene  distillate  as  medium  in  conjunction  with  a 
multiple  temperature  method,  namely  at  15°,  20°,  25°,  and  30  C,  so  that  the 
falling  drop  method  was  useful  for  analyzing  heavy  water  of  an  entire  concen¬ 
tration  range  with  an  accuracy  of  0.02-0.05  atom  percent  D. 

For  low  concentration  heavy  water  samples,  we  used  a  new  medium,  n-butyl 
benzoate,  and  studied  the  effect  of  the  residual  liquid  at  the  tip  of  the 
tube  during  "drop  making."  The  tube  was  treated  with  SiHCl3  so  that  no 
residual  liquid  stuck  to  the  tip.  We  also  studied  the  wall  effect  and  memory 
effect  and  confirmed  that  there  was  an  optimal  value  of  r/R,  the  ratio  of  drop 
radius  r  and  the  radius  of  the  falling. drop  tube  R,  when  the  drop  falling 
velocity  was  plotted  against  the  drop  volume.  The  drop  falling  velocity 
changed  very  little  at  the  vicinity  of  this  optimum.  The  optimum^determined 
for  n-butyl  benzoate  was  0.426  and  for  o-fluorotoluene  0.439.  »  »  There¬ 

fore,  it  is  much  more  convenient  to  use  the  popular  "microsampler  to  inject 
sample  drops  of  10.00±0.02  pi. 

Later,  a  mixture  of  a-ehloronaphthalene  and  n-butyl  benzoate  was  used  to 
analyze  samples  of  high  concentration  heavy  water  with  97  to  100  atom  percent 
D.  An  r/R  value  of  0.430  was  chosen.  The  reproducibility  of  the  drop  falling 
velocity  was  0.2  percent  (relative)  and  the  accuracy  of  the  determinations 
reached  ±0 . 04  atom  percent  D . 6 0 

3.5  The  Total  Isotope  Analysis  of  Water  by  the  Density  Method 

When  mass  spectrometer  is  not  available,  the  density  method  is  often  employed 
to  determine  the  total  content  of  D20  and  H2l80.  Then  D  and  0  are  normal¬ 
ized  and  one  of  them  measured.  The  difference  is  attributed  to  the  other. 

In  1936,  a  hydrogen  isotopes  normalization  method  using  H2S  was  studied  an 
a  C02  normalization  of  oxygen  isotopes  was  used  in  1959.  Later,  the  oxygen 
isotopes  normalization  was  done  over  Mn02  catalyst  at  600  C  using  compressed 
oxygen.  The  cold  trap  and  saturator  on  both  sides  of  the  horizontal  reaction 
chamber  were  of  the  same  shape  so  that  one  forward-reverse  operation  was. 
enough  to  normalize  the  unusually  high  concentration  of  oxygen  isotopes  in 
natural  water  to  within  the  error  range  of  measurement.  Since  the  1970  s, 
the  mass  spectrometry  has  been  used  for  determining  D,  0  and  0  separately. 

3.6  The  Mass  Spectrometric  Analysis  of  Hydrogen  Isotopes 

The  deuterium  concentration  of  a  heavy  water  product  was  examined  by  using  a 
MAT-CH5  mass  spectrometer.  In  actual  practice,  one  measured  t  e  ratio  o  pea 
height,  HD**"/DD+.  In  order  to  compensate  for  the  mass  discrimination  factor, 
a  reference  sample  of  145.7  ppmD  was  used  for  dilution.  We  compared  the  high 
temperature  platinum  wire  method  and  the  Pt02  catalysis  for  the  preparation  of 
gas  samples  and  found  the  latter  more  convenient.  When  calculating  molecular 
weights,  the  contents  of  170  and  180  have  to  be  taken  into  consideration^  The 
results  showed  the  purity  of  the  sample  was  99.939±0.003  atom  percent  D. 

This  sample  was  also  sent  to  the  labs  in  Japan  and  France  to  confirm  our 
results  and  the  data  were  in  agreement.  The  mass  spectrometry  has  been  used 
for  the  preparation  of  second  and  third  batches  of  standard  sample  series  so 
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that  a  standard  curve  can  be  obtained  for  production  departments5  to  test  their 
products  by  simple  infra-red  measurement. 

3.7  The  Reiteration  Method  in  the  Mass  Spectrometric  Analysis  of  Lithium 
Isotopes 

In  conjunction  with  the  research  on  Lithium-6,  we  studied  the  mass  spectro¬ 
metric  analysis  of  lithium  isotopes  in  1964. 50  A  MM-1305  mass  spectrometer 
was  used.  As  an  ion  source,  a  tantalum  strip,  after  being  applied  a  LiN(>3 
solution,  was  dried  on  hot  plate  to  form  a  solid  layer  of  600  pg/cm2.  When 
vaporized  at  900°C,  the  lithium  salt  emitted  6Li+  and  7Li+  ions  and  R7/6  was 
measured .  ' 

After  being  purified  by  crystallization,  carbonates  of  natural  and  enriched 
lithium  were  examined  for  impurities  and  found  to  be  mainly  potassium, 
sodium,  calcium,  and  iron.  An  intermediate  sample  was  obtained  by  dissolving 
and  thorough  mixing  of  these  two  samples  which  were  accurately  weighed.  By 
comparing  the  calculated  and  spectrometrically  determined  R7/6  values  of  the 
intermediate  sample,  a  first-approximation  correction  factor  . 

Kcalc./exp.  =  0.9944  was  obtained.  Through  four  reiterations,  the  values 
0.9926,  0.9920,  0.9918,  and  0.9918  were  obtained  successively.  By  using  this 
converged  K  value  to  do  corrections,  the  Li-6  concentrations  in  natural,  mixed, 
and  enriched  lithium  were  found  to  be  7.5,  49.2,  and  90.0  atom  percent  6Li, 
respectively.  The  standard  deviation  was  ±0.1  atom  percent  6Li.  The  enriched 
sample  Was  a  Russian— made  reagent.  The  6Li  concentration  of  the  mixed  sample 
was  also  checked  by  the  crystal  flotation  method  to  be  49.2±0.1  atom  percent 
6Li,  which  was  in  total  agreement. 

3.8  The  Elimination  of  Memory  Effect  in  the  Boron  Trifluoride  Method  for  Mass 
Spectrometric  Analysis 

In  connection  with  running  a  pilot  plant  for  enriching  B-10,  we  studied  the 
BF3  method  in  1965. 5 2  Because  the  memory  effect  in  this  method  is  mainly 
attributed  to  the  contamination  of  the  line,  the  BF3  samples  were  adsorbed  to 
pure  microcrystalline  CaF2  at  155°  to  form  complexes.  A  tantalum  or  nickel 
crucible  containing  10  mg  of  BF3*3CaF2  was  inserted  into  the  ion  source  inside 
a  Russian-made  MIT-1305  mass  spectrometer  and  heated  to  generate  BF 3  gas, 
which  was  then  bombarded  with  electrons  to  form  10BFf  and  11BFt  ions.  The 
R(11B/10B)  values  were  determined.  The  memory  effect  was  not  observed  when 
high  and  low  concentration  10B  samples  were  alternated  in  the  analysis 
sequence.  The  interference  of  other  ions  was  discussed  and  their  total  influ¬ 
ence  on  the  R  value  would  not  exceed  1  percent.  An  enriched  sample  was 
examined  by  this  method  and  found  to  have  82.2±0.1  atom  percent  10B,  which  was 
in  total  agreement  with  results  from  the  borax  method.  The  system  error 
correction  was  not  made  because  only  relative  measurements  were  made  for 
samples  taken  from  each  section  of  the  separation  column. 

3.9  The  Analysis  of  Reference  Samples  of  Nitrogen  Isotopes  by  Mass  Spectrometry 

Through  mass  discrimination  correction,  the  15N  abundances  of  three  ammonium 
sulfate  reagents  were  accurately  determined  by  using  a  MAT-251  mass 
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spectrometer.  Their  purities  were  also  determined  and  listed  in  Table  9. 
These  reagents  can  serve  as  reference  samples  for  the  preparation  of  inter¬ 
mediate  abundance  samples  useful  for  calibrating  mass  spectrometers  and  other 
instruments. 6 9  We  also  prepared  a  set  of  eight  low  15N  reference  samples  for 
agricultural  applications.  The  abundance  ranges  from  0.366  to  12.6  atom 
percent  1SN. 68 


Table  9. 

Reference 

Sample 

Purity, 

percent 

Atom 

percent  1SN 

I 

99.77 

0.36610.001 

II 

99.28 

98.3510.02 

III 

99.93 

99.8810.01 

3.10  The  Preparation  of  013  Reference  Samples  and  the  Mass  Spectrometric 
Analysis  of  Oxygen  Isotopes 

Like  the  nitrogen  isotope  project,  the  preparation  of  reference  samples  of 
carbon  isotopes  is  also  on  the  agenda  of  the  State  S&T  Commission.  A  99.95 
percent  pure  barium  carbonate  reagent  was  treated  with  100  percent  H3PO4  to 
get  C02  sample  gas  in  98.9±0.3  percent  conversion.  The  %  3/1 2  was  measured 
with  a  MAT-251  mass  spectrometer.  In  actual  practice,  three  peaks  at  44, 

45,  and  46  were  measured  (the  other  three  minor  peaks  were  ignored)  and 
corrected  by  the  170  and  180  content  calculations  to  get  Ri3/i2-  The  mass 
discrimination  factor  can  be  calculated  by  using  the  PDB  reference  samples 
for  comparison  with  the  <$13/12  and  Sia/i6  values  obtained  for  comptesse 
C02,  from  Which  the  true  value  of  R13/12  can  be  figured  out.  We  obtained 
the  result  of  1.095±0.0001  atom  percent  13C.  The  oxygen  isotope  contents 
were  also  determined,  i.e.  0.0375+0.0001  atom  percent  0  and  0.2032±0.0008 
atom  percent  0. 

Later ,  a  high  C-13  barium  carbonate  reagent  was  measured  and  its  isotope  abun¬ 
dance  of  88. 5410.03  atom  percent  13C  was  determined  after  the  correction  for 
mass  discrimination  factor  was  made.  The  0-18  in  the  reagent  was  depleted  to 
0.223  atom  percent  180.  We  also  found  a  reagent  containing  98.1110.03  atom 

percent  13C. 

IV.  The  Natural  Abundance  of  Hydrogen  and  Oxygen  Isotopes 

4.1  The  Distribution  of  Deuterium  in  the  Precipitation,  Surface  and  Under¬ 
ground  Waters 

At  the  end  of  1935,  four  snow  samples  were  collected  in  Northern  Europe6  and 
the  change  in  the  overall  content  of  D  and  180  measured  by  the  float  method. 

In  comparison  with  ocean  water,  the  density  difference  ranges  from  -3.7  to 
-4.9y  and  the  difference  is  -4.2y  when  compared  with  the  vapor  from  evaporation 

of  ocean  water  at  0°C.  ,  j  '■ 
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In  1959,  we  collected  all  kinds  of  water  samples  from  all  over  the  country 
including  waters  from  sea,  river,  lake,  underground  spring,  hot  spring, 
plateau  lake,  salt  lake,  and  oil  brine,  a  total  of  50  samples.  The  total 
content  of  D  and  180  was  analyzed  by  the  float  method  while  the  content  of 
D  and  180  of  a  few  samples  were  determined  by  the  C02  balance  method. 

The  more  obvious  results  were  that  ocean  water  was  heavier  than  surface 
water,  plateau  lake  water  heavier  than  plain  lake  water,  and  waters  of  salt 
lake  and  oil  brine  were  even  heavier.  This  can  generally  be  explained  by 
the  vapor  pressures  of  H20,  HDO ,  and  H2180.37 

Between  1956  and  1983  we  collected  six  batches  of  samples  from  the  Chang 
Jiang  (Wuhan,  Nanjing,  and  Shanghai),  the  Zhu  Jiang  (Guangzhou),  and  the 
Qian  Tang  Jiang  (Hangzhou) .  The  waters  were  collected  for  12  consecutive 
months.  Reduction  with  magnesium  powder  produced  hydrogen  gases  for  precise 
chromatography.  The  results  showed  that  deuterium  content  was  highest  in 
May  and  lowest  in  September,  the  difference  being  about  4-5  ppmD.  This  was 
consistent  with  previous  results.  The  difference  is  due  to  the  fact  that  ' 
river  water  comes  from  precipitations,  which  are  gathered  from  vapors 
evaporized  at  different  temperatures. .  At  higher  temperature,  the  ratio  of 
vapor  pressure  HD0/H20  (less  than  1  itself)  deviates  less  from  unity. 

Between  1974  and  1983  we  collected  the  deep  well  waters  of  the  Beijing  Univer¬ 
sity,  including  the  12  months  samples  of  1978  and  1981-82.  The  deuterium 
contents  were  measured  as  145.7±0.2  ppmD,  an  indication  that  it  was  constant 
within  the  error  range  of  the  determination  in  a  9-year  span.70 

In  the  1970's,  we  measured  the  hydrogen  and  oxygen  isotopes  distribution  of 
snow  and  ice  samples  in  the  high  altitude  area  of  Mt.  Qomolangma. 4 6-4 8  The 
density  method  for  the  analysis  of  total  isotopes  and  the  mass  spectrometry 
were  used  and  some  data  are  listed  in  Table  10.  These  samples  were  collected 
between  21-27  April  1975. 


Table  10. 


Attitude , 
m 

Precipitation 

<$D  (%0) 

<5 1 70  (%c) 

<5180  (%0) 

7450 

Granular  snow 

-  62.0 

-  5.25 

-  9.43 

Ice 

-  22.3 

-  1.86 

-  4.31 

Granular  snow 

-  82.2 

-  6.23 

-12.57 

Surface  snow 

-136.2 

-  9.17 

-18.79 

8848 

New  snow 

-208.8 

-13.82 

-26.99 

4.2  The  Regularity  in  the  Distribution  of  Hydrogen  and  Oxygen  Isotopes  on 
Earth 

In  1936,  we  proposed  that  oxygen  and  hydrogen  isotopes  recycled  in  nature  in 
the  form  of  water  based  on  the  fact. that  the  vapor  pressures  of  HDO  and  H2180 
are  lower  than  that  of  H20.  Water  vaporizes  and  rises  to  form  cloud,  which 
in  turn  condenses  or  freezes  and  descends  in  the  form  of  rain  and  snow.  Part 
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of  the  vapors  from  ocean  water  precipitate  on  land  and  surface  waters  flow 
one-way  into  the  ocean,  hence  forms  a  cycle  and  D  and  0  are  enriched^ in 
ocean  waters  (on  average,  surface  water  contains  148  ppmD  and  1970  ppm  0 
whereas  ocean  water  156  ppmD  and  2000  ppmle0).  At  the  same  time,  land  animals 
also  enrich  heavy  isotopes  from  the  water  they  drink  through  surface  evapora¬ 
tion.  Plants,  due  to  their  larger  volume  of  evaporation,  have  greater  enrich¬ 
ment  (the  crystallization  water  of  minerals  also  has  enrichment  effect  compar¬ 
ing  with  ocean  water).  In  conclusion,  the  evaporation  formulas  of  HDO  and 
H2180  can  be  used  to  roughly  explain  the  distribution  of  D  and  180  on  earth, 
although  there  are  other  secondary  factors . 7 »8 »9 

4.3  The  Effect  of  the  Fluctuation  of  Hydrogen  and  Oxygen  Distribution  on  the 
Units  of  Volume  and  Temperature 

In  1936,  we  discussed  the  problem  of  such  metric  standards  as  volume  and 
temperature10  and  pointed  out  that,  in  measuring  1  liter  standard  volume  by 
water,  the  source  of  the  regular  water  and  its  purification  method  should  be 
specified.  It  is  still  impractical  to  use  the  ideal  H2  0  as  the  standard 
material.  The  fluctuation  of  hydrogen  and  oxygen  isotopes  has  smaller  effect 
on  temperature  measurements.  With  the  new  international  definition  of 
ll=ldm3 ,  this  problem  is  solved. 

Looking  back  on  my  50  years  of  experience  in  scientific  research,  I  spent  the 
early  years  on  basic  research,  the  results  of  which  have  been  applied  to 
separation  and  analysis  methodology.  After  the  Liberation,  in  adjusting  to 
the  demands  of  our  socialist  construction,  I  turned  to  applied  research, 
which  was  focused  on  the  study  of  chemistry  and  analysis  of  stable  isotopes 
of  light  elements  with  an  aim  for  their  production  and  development.  Half 
century  of  efforts  is  just  like  a  drop  in  the  bucket.  But  its  path  truly 
reflects  the  direction  the  surging  torrent  of  new  China's  S&T  development  is 
heading. 


REFERENCES 

1.  Riesenfeld,  E.H. ,  Zhang  Qinglian,  Z.  PHYSIK.  CHEM. ,  B28,  408-409  (1935). 

2.  Ibid.,  B30,  61-68  (1935). 

3.  Ibid.,  B33,  120-126  (1936). 

4.  Ibid.,  127-132. 

5.  Riesenfeld,  E.H.,  Zhang  Qinglian,  BER.  DTSCH.  CHEM.  GES. ,  69,  1302-1305 
(1936). 

6.  Ibid.,  1305-1307. 

7.  Ibid.,  1308-1310. 

8.  Riesenfeld,  E.H. ,  Zhang  Qinglian,  NATURWISS. ,  24,  616-618  (1936). 


75 


9. 


Riesenf eld ,  E.H. ,  Zhang  Qinglian,  EMSCHAU  WISS.  TECHNIK,  40,  621-623 


10.  Riesenfeld ,  Zhang  Qinglian,  PHYSIK,  Z.,  37,  690-692  (1936). 

11.  Zhang  Qinglian,  "The  Isotope  Chemistry  of  Water,"  Diss.  Univ.  Berlin, 

50  pages,  Neuenhahn,  Jena,  1936. 

12.  Zhang  Qinglian,  Zhu  Jinchang,  Z.  PHYSIK.  CHEM. ,  A184,  411-415  (1939). 

13.  Zhang  Qinglian,  Wei  Yuzhi,  ZH0NGGU0  HUAXUEHUI  HUIZHI  [JOURNAL  OF 
THE  CHINESE  CHEMICAL  SOCIETY],  7,  138-143  (1940). 

14.  Zhang  Qinglian,  Qian  Renyuan,  Ibid.,  8,  74-75  (1941). 

15.  Zhang  Qinglian,  Qian  Renyuan,  J.  AM.  CHEM.  S0C.,  63,  1709-1711  (1941). 

16.  Zhang  Qinglian,  Wei  Yuzhi,  KEXUE  JILU  [SCIENTIFIC  RECORDS],  1,  132-134 
(1942). 

17.  Zhang  Qinglian,  SCIENCE,  100,  29-30  (1944). 

18.  Zhang  Qinglian,  J.  AM.  CHEM.  SOC.,  66,  1940-1941  (1944). 

19.  Zhang  Qinglian,  Zhu  Jinchang,  QINGHUA  DAXUE  LIKE  BAOGAO  [PHYSICAL 
SCIENCES  REPORT,  QINGHUA  UNIVERSITY],  A4,  361-362  (1947). 

20.  Zhang  Qinglian,  Wei  Yuzhi,  Ibid.,  A4,  451-452  (1947). 

21.  Zhang  Qinglian,  ZH0NGGU0  KEXUE  JISHU  [CHINESE  SCIENCE  AND  TECHNOLOGY],  1, 
109-110  (1948) . 

22.  Zhang  Qinglian,  "Studies  on  Heavy  Water,"  the  50th  Anniversary  Commemora¬ 
tive  Publication  of  Beijing  University,  The  School  of  Physical  Sciences, 
129-142  (1948). 

23.  Zhang  Qinglian,  Xie  Yuyuan,  QINGHUA  DAXUE  LIKE  BAOGAO  [PHYSICAL  SCIENCES 
REPORT,  QINGHUA  UNIVERSITY],  A5 ,  252-259  (1948). 

24.  Zhang  Qinglian,  Xie  Yuyuan,  ZH0NGGU0  HUAXUEHUI  HUIZHI  [JOURNAL  OF  THE 
CHINESE  CHEMICAL  SOCIETY],  16,  10-13  (1949). 

25.  Ibid.,  65-71. 

26.  Zhang  Qinglian ,  Dong  Luhe,  NATURE,  163,  737  (1949). 

27.  Zhang  Qinglian,  Dong  Luhe,  ZHONGGUO  WULI  XUEBAO  [ACTA  PHYSICA  SINICA] , 

7,  230-240  (1949). 

28.  Zhang  Qinglian,  Dong  Luhe,  ZHONGGUO  HUAXUEHUI  HUIZHI  [JOURNAL  OF  THE 
CHINESE  CHEMICAL  SOCIETY],  16,  1-9  (1949). 


76 


29.  Zhang  Qinglian,  Wang  Jinhua,  Ibid.,  59-64. 

30.  Zhang  Qinglian,  Zeng  Yinglun,  KEXUE  JILU  [SCIENTIFIC  RECORDS],  3,  101-105 
(1950). 

31.  Zhang  Qinglian,  ZHONGGUO  KEXUE  [SCIENTIA  SINICA] ,  2,  185—190  (1951). 

32.  Zhang  Qinglian,  Zhang  Rongsen,  HUAXUE  XUEBAO  [ACTA  CHIMICA  SINICA],  21, 
320-327  (1955);  ZHONGGUO  KEXUE  [SCIENTIA  SINICA],  4,  555-561  (1955). 

33.  Zhang  Qinglian,  Gao  Zi,  Huang  Chunhui,  Ibid.,  22,  234-240  (1956). 

34.  Zhang  Qinglian,  Wang  Kechang,  Ibid.,  22,  414-416  (1956). 

35.  Zhang  Qinglian,  Cheng  Zhixiang,  Xu  Huibi,  KEXUE  JILU  [SCIENTIFIC  RECORDS, 
NEW  EDITION],  1,  383-384  (1957);  English  Edition,  411-413  (1957). 

36.  Zhang  Qinglian,  Cheng  Zhixiang,  Xu  Huibi,  "Production  of  Isotopes,  etc." 
IZV  AN  SSSR,  59-66  (1958),  paper,  Conference  on  Isotope  Applications, 
Moscow,  April  1957. 

37.  Zhang  Qinglian,  Yin  Zhao,  Zhang  Xianlu,  KEXUE  JILU  [SCIENTIFIC  RECORDS, 
NEW  EDITION],  3,  402-404  (1959);  English  Edition,  495-498  (1959). 

38.  Zhang  Qinglian,  "Studies  on  the  Heavy  Water  Analysis,"  Higher  Education 
Publication,  p  22  (1959). 

39.  Huang  Ailun,  Guo  Dezhu,  Jiao  Rangjie,  Cheng  Zhixiang,  Zhang  Qinglian, 
KEXUE  JILU  [SCIENTIFIC  RECORDS,  NEW  EDITION] ,  3,  189-191  (1959);  English 
Edition,  239-245  (1959). 

40.  Zhang  Qinglian,  KERNENERGIE ,  3,  892-896  (1960),  paper,  Conference  on 
Isotope  Applications,  Leipzig,  December  1959. 

41.  Zhang  Qinglian ,  "The  Physical  Properties  of  Heavy  Water  and  Heavy 
Hydrogen,"  document,  820  Professional  Meeting,  p  34,  September  1960. 

42.  Zhang  Qinglian,  "Methods  for  Heavy  Water  Analysis,"  Ibid.,  p  29. 

43.  Jin  Deqiu,  Zhang  Qinglian,  ZHONGGUO  KEXUE  [SCIENTIA  SINICA],  12, 

147-148  (1963). 

44.  Jin  Deqiu,  Zhang  Qinglian,  HUAXUE  XUEBAO  [ACTA  CHIMICA  SINICA],  29, 
266-276  (1963). 

45.  Zhang  Qinglian,  Gu  Zhennan,  ZHONGGUO  KEXUE  [SCIENTIA  SINICA],  13, 

688-689  (1964). 

46.  Zhang  Shen,  Yu  Weixin,  Zhang  Qinglian,  Huang  Chunhui,  Zhang  Rongsen, 

Gu  Zhennan,  Wei  Qingong,  Xie  Zichu,  Zeng  Qunshu,  ZHONGGUO  KEXUE  [SCIENTIA 
SINICA],  16,  430-433  (1973);  English  Edition,  560-564  (1973). 


77 


47.  Zhang  Shen,  Zhang  Qinglian,  Xie  Zichu,  "Report  of  Scientific  Survey  in 
Mt.  Qomolangma  Area  (1966-1968),  Modern  Glaciers  and  Geomorphology," 
Science  Publication,  71-76  (1975). 

48.  Zhang  Qinglian, >  KEXUE  TONGBAO  [MONTHLY  JOURNAL  OF  SCIENCE],  English 
Edition,  25,  208-212  (1980),  paper,  27th  IUPAC  Congress,  Helsinki, 
August  1979. 

49.  Zhang  Qinglian,  Huang  Hantong,  SEPAR.  SCI.  TECHN.,  15,  397-403  (1980), 
abstracts,  p  3,  Separation  Science  and  Technology  Symposium,  Gatlinburg, 
October-November  1979. 


50.  Zhang  Qinglian,  Wang  Guanghui,  Zhuan  Yahui,  KEXUE  TONGBAO  [MONTHLY 
JOURNAL  OF  SCIENCE],  English  Edition,  26,  897-900  (1981),  paper,  Gordon 
Conference  on  Isotope  Chemistry  and  Physics,  Ventura,  February  1981. 

51.  Zhang  Qinglian ,  Shen  Dixin,  Ibid.,  26,  995-1000  (1981). 

52.  Wang  Guanghui,  Zhang  Qinglian,  Ibid.,  27,  986-988  (1982);  English 
Edition,  968-972  (1982),  12th  British  Conference  of  Mass  Spectroscopy, 
Cambridge,  September  1981. 

53.  Zhou  Xihuang ,  Xu  Dagang,  Xu  Wenrui,  Zhang  Qinglian,  HUAXUE  TONGBAO 
[CHEMISTRY],  (9),  568-569  (1982). 

54.  Zhou  Xihuang,  Li  Ping,  Yang  Guomin,  Zhang  Qinglian,  Ibid.,  (6),  17-18 


55.  Yu  Jianguo,  Ni  Baoling,  Zhang  Qinglian,  KEXUE  TONGBAO  [MONTHLY  JOURNAL 
OF  SCIENCE],  28,  405-407  (1983);  English  Edition,  29,  776-779  (1984). 

56.  Nan  Chunbo,  Qian  Qiuyu,  Zhang  Qinglian,  ZH0NGGU0  KEXUE  [SCIENTIA  SINICA] 
English  Edition,  B26 ,  449-453  (1983). 


57. 


Huang  Chunhui,  Zhang  Qinglian,  HUAXUE  TONGBAO  [CHEMISTRY],  (8).  11-12 
(1983). 


58.  Ibid.,  12-13. 

59.  Wang  Shouguo,  Zhang  Qinglian,  Ibid.,  (9),  17-18  (1983). 

60.  Ibid. ,18-19. 

61.  Zhang  Qinglian,  "Stable  Isotopes  Research  at  the  Beijing  University," 
paper,  Saclay  Atomic  Nucleus  Research  Center  Workshop,  August  1983. 

62.  He  Shaoqi,  Wu  Jinguang,  Zhang  Qinglian,  REV.  CHIM.  MINTR. ,  20,  737-744 


63.  He  Shaoqi,  Wu  Jinguang,  Zhang  Qinglian,  HUAXUE  XUEBAO  [ACTA  CHIMICA 
SINICA],  42,  1183-1187  (1984). 


78 


64.  Zhou  Xihuang ,  He  Xihua,  Zhang  Qinglian,  HUAXUE  TONGBAO  [CHEMISTRY],  (1) , 
14-15  (1984). 

65.  Qian  Qiuyu,  Wang  Mei,  Li  Hong,  Zhang  Qinglian,  KEXUE  TONGBAO  [MONTHLY 
JOURNAL  OF  SCIENCE],  29,  402-403  (1984);  English  Edition,  897-899  (1985) 

66.  Qian  Qiuyu,  Qian  Jiuxin,  Wang  Baiyang,  Zhang  Qinglian,  Ibid.,  785-788. 

67.  Qian  Giuyu,  Zhang  Qinglian,  Ibid.,  1058-1061. 

68.  Yang  Yuemei,  Jin  Deqiu,  Zhang  Qinglian,  Ibid.,  1504-1506. 

69.  Yang  Yuemei:,  Jin  Deqiu,  Zhang  Qinglian,  "Proceeding  of  Sino-Japanese 
Mass  Spectroscopy  Symposium,"  English  Edition,  Tokyo,  1985. 

70.  Zhou  Xihuang,  Jiang  Shihe,  Zhang  Qinglian,  HUAXUE  TONGBAO  [CHEMISTRY], 
(1),  13-14  (1985). 

71.  Li  Wenjun,  Qian  Qiuyu,  Zhang  Qinglian,  Ibid.,  (7),  23  (1985). 

72.  Qian  Qiuyu,  Qian  Jiuxin,  Ye  Yijun,  Zhang  Qinglian,  Ibid.,  (6),  23  (1985) 

12922/9365 
CSO:  4008/1086 


79 


Acoustics 

JPRS*CST*85*040 
20  November  1985 

AUTHOR:  WANG  Shuozhong  [3769  2592  0022] 

ORG:  Institute  of  Acoustics,  Chinese  Academy  of  Sciences 

TITLE:  "An  Application  of  Diffraction  Theory  to  Echo  Formation" 

SOURCE:  Beijing  SHENGXUE  XUEBAO  [ACTA  ACUSTICA]  in  Chinese  Vol  10  No  5 
Sep  85  pp  302-318 

TEXT  OF  ENGLISH  ABSTRACT:  Diffraction  theory  is  applied  to  treat  echo- 
formation  problems  in  this  paper.  This  overcomes  some  drawbacks  of  the 
previously  accepted  image  pulse  theory.  The  new  approach  provides  a 
comprehensive  physical  interpretation  for  the  mechanism  of  echo  formation. 

It  predicts  different  echo  behavior  for  objects  made  of  different  materials, 
and  gives  correct  results  for  problems  of  echo  formation  by  smoothly  curved 
surfaces.  In  these  aspects,  the  previous  theory  breaks  down.  In  principle, 
the  diffraction  approach  can  tackle  any  echo-formation  problem  so  long  as 
there  exists  a  solution  to  the  relevant  configuration  of  diffraction.  The 
method  is  strongly  supported  by  experimental  results. 
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Sep  85  pp  319-323 

TEXT  OF  ENGLISH  ABSTRACT:  The  formula  for  sound  velocity  in  a  mixture  of 
weak  interactives  is  used  to  calculate  the  sound  velocities  in  seven  ternary 
systems.  The  results  show  that  the  agreement  between  the  calculated  and 
experimental  values  of  sound  velocity  is  quite  good.  The  results  illustrated 
in  this  article  show  that  the  large  negative  value  of  excess  isentropic 
compressibility  does  not  always  indicate  a  large  positive  value  of  (9VE/9p)s. 
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TITLE:  "The  Experimental  Investigation  of  Karman  Vortex  Street  in  the 

Supersonic  Wake" 

SOURCE:  Mianyang  KONGQIDONGLIXUE  XUEBA0  [ACTA  AERODYNAMICA  SINICA] 
in  Chinese  No  3,  1985  pp  1-8 

TEXT  OF  ENGLISH  ABSTRACT:  A  Karman  vortex  street  was  observed  in  the  wake  of 
a  rough  flat  plate  at  M  =  2.2.  The  supersonic  vortex  flow  was  investigated 
by  optical  methods.  Flow  velocity  and  velocity  fluctuations  were  recorded 
with  a  laser  doppler  anemometer  (LDA) .  The  flow  was  visualized  by  schlieren 
and  interferometer  methods,  and  the  density  distribution  in  the  wake  was 
measured  with  a  Mach-Zehnder  interferometer  in  connection  with  a  high  speed 
camera.  The  frequency  of  vortex  shedding  was  determined  with  a  laser 
interferometer.  ■ 


The  width  of  the  wake  is  larger  and  base  pressure  of  the  plate  is  lower  when 
compared  with  a  fully  turbulent  wake  behind  a  smooth  flat  plate  under  the  same 
flow  conditions.  The  ratio  h/1  of  the  vortex  street  is  about  0.37,  i.e., 
higher  than  that  in  the  incompressible  flow  behind  a  cylinder.  The  Strouhal 
number  of  the  vortex  street  shedding  is  about  0.22.  The  influence  of  the 
boundary  conditions  on  the  vortex  street  and  the  interaction  of  the  vortex 
street  with  an  oblique  shock  are  discussed  in  the  paper. 
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TEXT  OF  ENGLISH  ABSTRACT:  Based  on  the  separated  vortex-field  of  the  wing 
and  vortex-lifting  theory,  the  theoretical  formulas  of  the  supersonic  wing's 
rolling  vortex-damping  due  to  the  separated  vortex  of  the  low-aspect-ratio 
sharp-edged  wings  under  the  existence  of  the  roll  and  attack  angle  are 
derived  using  supersonic  linear  theory.  Two  terms  are  contained  in  the  roll 
vortex  damping — one  is  proportional  to  the  product  of  the  attack  angle  and 
roll  rate,  and  the  other  is  proportional  to  the  square  of  the  roll  rate. 
Then,  the  new  theoretical  calculation  model  of  the  roll  damping  is  given. 

From  the  theoretical  results  obtained,  some  figures  of  the  rolling  damping 
factor  with  various  combinations  of  parameters  of  wing  platforms  are  also 
given.  In  addition,  some  comparisons  between  this  method  and  the  prevalent 
potential  theory  are  presented. 
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Separating  Vortex  for  Slender  Delta  Wings" 

SOURCE:  Mianyang  KONGQIDONGLIXUE  XUEBAO  [ACTA  AERODYNAMICA  SINICA] 
in  Chinese  No  3,  1985  pp  63-71 

TEXT  OF  ENGLISH  ABSTRACT:  A  double  layer  mode  was  established,  which  is  used 
here  to  describe  the  separating  vortex  flow  over  the  delta  wing  leading-edge. 
The  inner  layer  of  this  double  layer  mode  is  the  vortex  core  and  the  outer 
layer  is  the  irrotational  potential  flow.  These  provide  the  initial  and  outer 
edge  conditions  which  are  necessary  for  solving  vortex  core  development.  In 
incompressible  flow  this  method  can  not  only  compute  the  initial  breakdown 
angles  of  attack,  but  can  also  compute  the  breakdown  points  along  the  line  of 
the  root  chord  when  the  angles  of  attack  are  greater  than  the  critical  ones. 
The  theoretical  results  agree  well  with  the  experimental  measures  and  show 
that  the  breakdown  points  are  insensitive  to  the  Reynolds  number. 
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TITLE:  "On  Rolling  Resonance  and  Small  Asymmetrical  Aerodynamic  Force  of 

Re-entry  Body" 

SOURCE:  Mianyang  KONGQIDONGLIXUE  XUEBAO  [ACTA  AERODYNAMICA  SINICA] 
in  Chinese  No  3,  1985  pp  72-81 

TEXT  OF  ENGLISH  ABSTRACT:  In  this  paper,  the  problems  in  flight  mechanics 
and  aerodynamics  of  rolling  resonance  of  the  re-entry  body  are  discussed;  the 
mechanism  of  the  divergence  of  the  angle  of  attack  of  the  re-entry  body  is 
also  analyzed.  It  is  shown  that  rolling  resonance  is  a  primary  hazard  of  the 
re-entry  body  in  the  divergence  of  the  angle  of  attack.  Research  results  of 
the  rolling  resonance'  and  small  asymmetrical  aerodynamic  force  as  well  as 
existential  problems  relating  to  the  re-entry  body  are  emphatically  described. 
It  is  shown  that  the  determination  of  the  random  ablative  configuration  and  the 
corresponding  mathematical  expectation  and  the  variance  of  small  asymmetrical 
aerodynamic  coefficients  is  the  key  problem  in  predicting  the  rolling 
resonance  probability  of  the  re-entry  body.  Flight  testing  is  the  primary 
method  for  solving  this  difficult  problem. 
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SCOTT,  S.K. 

ORG:  SCOTT  of  the  School  of  Chemistry,  University  of  Leeds,  United  Kingdom 

TITLE:  "Thermal  Explosion  of  Dispersed  Media:  Criticality  for  Discrete 

Reactive  Particles  in  an  Inert  Matrix" 

SOURCE:  Beijing  BINGGONG  XUEBAO  [ACTA  ARMAMENT ARI I ]  in  Chinese  No  3,  Aug  85 

pp  1-10 

TEXT  OF  ENGLISH  ABSTRACT:  In  the  present  paper,  theoretical  interpretations 
are  developed  which  explore  the  criteria  for  safe  storage  of  explosive  masses 
in  the  situation  in  which  a  number  of  centers  of  reaction  are  embedded  in  an 
inert,  but  not  perfectly  conducting,  matrix.  The  simplest  case  is  that  in 
which  the  temperature  of  the  matrix  is  uniform,  but  temperature  distributions 
within  the  reactive  zones  are  also  considered.  Finally  different  applications 
of  model  propellant  systems  are  discussed. 
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AUTHOR:  WU  Xiong  [0702  7160] 

ORG:  None 

TITLE:  "Detonation  Properties  of  Condensed  Explosives  Compacted  with  the 

VLW  Equation  of  State" 

SOURCE:  Beijing  BINGGONG  XUEBAO  [ACTA  ARMAMENTARII]  in  Chinese  No  3,  Aug  85 

pp  11-20 

TEXT  OF  ENGLISH  ABSTRACT:  A  simplified  virial  equation  of  state  based  on 
Lennard-Jones  6-12  potential,  called  VLW  EOS,  is  proposed  for  detonation 
products.  The  detonation  properties  of  explosives  composed  of  CHNO,  CNO, 

HNO,  NO,  HN,  CHNOF  and  CHNF  have  been  calculated  with  this  equation,  and  the 
results  are  compared  with  those  of  BKW  EOS  and  LJD  EOS.  The  principal 
feature  of  this  approach  is  that  the  third  and  fourth  virial  coefficients 
can  be  conveniently  obtained  from  the  well  known  second  one.  It  is  shown  that 
the  agreement  between  the  predicted  detonation  parameters  and  the  experimental 
data  is  reasonably  satisfactory. 
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ORG:  None 

TITLE:  "Design  of  Linear  Prismatic  Scanner  with  High  Accuracy" 

SOURCE:  Beijing  BINGGONG  XUEBAO  [ACTA  ARMAMENTARI I ]  in  Chinese  No  3,  Aug  85 
pp  32-38 

TEXT  OF  ENGLISH  ABSTRACT:  This  paper  presents  the  design  of  a  high  accuracy 
linear  scanner  used  in  laser  beam  scanning  systems  for  automatic  inspection 
of  object  size  by  means  of  a  rectangular  prism  and  correct  selection  of  the 
rotation  center.  The  transmitted  laser  beam  scans  in  the  direction  precisely 
parallel  to  the  outlet  direction.  For  prismatic  glass  with  a  refraction  index 
of  about  1.9  and  a  measurement  region  of  up  to  43.74  mm,  the  maximum  linear 
error  is  less  than  ±20  nm.  When  the  measurement  region  is  limited  to  within 
39  mm,  the  maximum  linear  error  is  less  than  ±6  nm. 
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WANG  Cenghui  [3769  2582  6540] 

TIAN  Deyu  [3944  1795  0151] 

LI  Guozhen  [2621  0948  6966] 

XU  Liqing  [1776  7787  3237] 

ORG:  None 

TITLE:  "Synthesis  and  Properties  of  2,3,5,6-Tetranitrato-l,4-Dinitropiper- 

azine" 

SOURCE:  Beijing  BINGGONG  XUEBAO  [ACTA  ARMAMENTARII]  in  Chinese  No  3,  Aug  85 
pp  60-63 

TEXT  OF  ENGLISH  ABSTRACT:  A  new  energetics  compound,  2,3,5,6-tetranitrato- 
1,4-dinitropiperazine  (TN),  is  synthesized  and  its  properties  are  determined. 
The  method  for  preparing  intermediate  disodium  2,3,5,6-tetrahydroxy- 
piperazine-l,4-disulphonate  (SP)  and  2,5-diacet6xy-3,6-dinitrato-i,4-dinitro- 
piperazine  (DN)  for  TN  has  been  improved. 
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TITLE:  "Design  and  Microprocessor- Implementation  of  QMF  Used  for  SBC" 

SOURCE:  Shenyang  DONGBEI  GONGXUEYUAN  XUEBAO  [JOURNAL  OF  NORTHEAST  INSTITUTE 

OF  TECHNOLOGY]  in  Chinese  No  3,  Sep  85  pp  107-112 

ABSTRACT:  Using  the  prevailing  and  low-priced  microprocessors,  a  quadrature 

mirror  filter  (QMF)  system  is  designed  and  implemented  for  subband  coder  (SBC). 
Methods  are  proposed  to  improve  the  low  precision  and  processing  speed  of  the 
microprocessor,  such  as  the  distributed  algorithm,  additional  hardware  logic 
and  multimicroprocessor  structure.  A  sampling  frequency  up  to  9000  Hz  is 
attained  by  the  use  of  a  32  tap  FIR  filter  and  a  Z80A  CPU  if  the  resolution 
of  intermediate  stored  value  is  16-bit  and  that  of  input  signal  12-bit.  It 
is  sufficient  to  meet  the  requirement  of  real-time  speech  signal  processing. 
(Paper  received  24  October  1984) 
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TITLE:  "Eiesign  and  Microprocessor- Implementation  of  ADPCM  Encoder  for  SBC" 

SOURCE:  Shenyang  DONGBEI  GONGXUEYUAN  XUEBAO  [JOURNAL  OF  NORTHEAST  INSTITUTE 

FOR  TECHNOLOGY]  in  Chinese  No  3,  Sep  85  pp  113-119 

ABSTRACT:  Designs  and  implements  an  adaptive  differential  pulse-code¬ 

modulation  (ADPCM)  system  for  encoding  and  decoding  in  the  sub-band  (SBC) 
system  of  speech  based  on  a  Z-80  microprocessor.  At  the  transmitting  end,  the 
16-bit  binary  digital  signals  are  compressed  as  3  or  4-bit  codewords  by  an 
encoder ,  and  then  at  the  receiving  end ,  a  decoder  is  used  to  recover  these 
codewords.  With  optimal  parameters  selected,  the  system  offers  a  dynamic 
range  of  40  db  and  the  process  signal- to-noise  ratios  attainable  are  25  and 
33  db,  respectively,  for  the  3  and  4-bit  codings.  Implementing  with  no  more 
than  20  IC  chips  and  taking  a  Z-80  CPU  which  operates  at  4  MHz  clock 
frequency  as  the  ALU,  the  encoder  can  process  each  sample  value  in  a  duration 
about  from  3000  to  3200  CPU  clock  cycles  with  an  acceptable  input  signal 
frequency  greater  than  1200  Hz.  The  system  is  thus  available  to  the  real-time 
coding  of  speech  at  8  sub-bands  (with  sampling  frequency  8000  Hz).  With  the 
addition  of  multiplier  hardware,  the  processing  speed  of  encoder  can  be 
increased  in  so  much  that  the  acceptable  input  signal  sampling  frequency  is 
greater  than  4000  Hz  and  the  system  can  thus  be  used  for  the  real-time  coding 
of  speech  at  two  sub-bands.  Discusses  in  detail  the  algorithm  design,  step- 
size  parameters  and  IC  implementation  and  analyzes  the  microcomputer  simulation. 
(Paper  received  22  October  1984) 
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TITLE:  "The  Application  of  a  Microprocessor  in  the  Steam  Flow  Measuring 
System" 

SOURCE:  Shanghai  ZIDONGHUA  YIBIAO  [PROCESS  AUTOMATION  INSTRUMENTATION]  in 
Chinese  Vol  6,  No  9,  20  Sep  85  pp  28-32,  46 

TEXT  OF  ENGLISH  ABSTRACT:  In  the  process  of  steam  flow  metering,  more  and  more 
attention  has  been  devoted  to  the  compensation  of  temperature  and  pressure.  A 
metering  system  consisting  of  TP801  single-board  processor  has  been  described 
and  its  functions  and  working  principle  are  discussed.  The  strategy  of  the 
recommended  system  has  been  field  tested  and  proven  to  be  practical  and 
feasible. 
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SOURCE:  Shanghai  ZIDONGHUA  YIB1A0  [PROCESS-  AUTOMATION  INSTRUMENTATION]  in 
Chinese  Vol  6,  No  9,  20  Sep  85  pp  33-35,  27,  32,  46 

f:  i  ^ 

TEXT  OF  ENGLISH  ABSTRACT:  An  operator  panel  is  the  window  fob  the  operator  to 
monitor  and  run  the  process  as  well  as  the'  process  of  man-machine  interface 
between  the  operator  and  the  computer.  A  functional  key  panel  for  process 
control  system  is  presented  with  detailed  descriptions  of  its  functions  and 
key-controlled  service  program. 
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SOURCE:  Shanghai  ZIDONGHUA  YIBIAO  [PROCESS  AUTOMATION  INSTRUMENTATION]  in 
Chinese  Vol  6,  No  9,  20  Sep  85  pp  11-14,  18,  46 

TEXT  OF  ENGLISH  ABSTRACT:  A  miniature  eddy-current  transducer  based  on  high 
frequency  eddy-current  effect  is  described.  The  diameter  of  the  probe  is  only 
2.2-7  mm,  which  makes  it  especially  suitable  for  measuring  the  displacement 
and  vibration  of  minute  parts.  The  construction,  working  principle  and  main 
technical  data  of  the  transducer  are  discussed  in  full  detail.  Brief  examples 
of  some  applications  are  given. 
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TITLE:  "Radon  and  Dust  Reduction  and  Environmental  Protection  at  Uranium 

Open-pit  Mine" 

SOURCE:  Beijing  YOU  KUANG  YE  [URANIUM  MINING  AND  METALLURGY]  in  Chinese  Vol  4, 
No  3>  Aug  85  pp  6-12 

TEXT  OF  ENGLISH  ABSTRACT:  At  a  certain  uranium  open-pit  mine,  the  sources  of 
dust  and  radon  have  been  investigated  since  its  establishment.  In  places 
where  dust  and  radon  are  found*  a  series  of  safety  measures  and  safeguards  has 
been  adopted  to  reduce  dust  and  radon  emission  and  ensure  workers'  health  in 
uranium  mining.  Mine  water  and  waste  have  been  treated  appropriately  for  the 
goal  of  environmental  protection.  It  follows  that  the  surroundings  and 
natural  environment  of  the  uranium  mine  have  not  yet  been  appreciably  polluted 
by  radioactive  fallout. 
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TITLE:  "Mineral  Processing  for  Some  Uranium  Ore  of  Sedimentary  Deposit 

Containing  High  Carbonates  and  Sulphides" 

SOURCE:  Beijing  YOU  KUANG  YE  [URANIUM  MINING  AND  METALLURGY]  in  Chinese  Vol  4, 
No  3,  Aug  85  pp  13-18 

TEXT  OF  ENGLISH  ABSTRACT:  Either  carbonate  leaching  or  acid  leaching  of  the 
uranium  ore  of  some  sedimentary  deposit  containing  high  carbonates  and 
sulphides  could  not  give  satisfying  technical  and  economical  performances. 
Thus  a  flotation-leaching  process  was  proposed  and  described  therewith  in  this 
paper. 

First  of  all,  the  selective  flotation  of  the  pyrite  in  alkaline  medium  by 
isopropyl  xanthate  was  conducted.  The  sulphide  concentrate  containing  39.1 
percent  S  with  the  S  recovery  of  80.98  percent  was  obtained.  Subsequently, 
the  flotation  of  the  carbonaceous  matter  by  pine  oil-kerosene  was  carried  out. 
Finally,  the  thickened  reground  tailings  were  further  treated  by  flotation 
with  oxidized  paraffin  soap  as  the  collecting  reagent  for  carbonates.  The 
floated  carbonates  containing  0.013“0.015  U  could  be  rejected.  The  uranium 
containing  sulphide  concentrate  and  the  mixed  concentrates  were  separately 
treated  by  acid  leaching  with  the  consumption  of  66  kg  sulfuric  acid  per  ton 
of  the  ore  subject  to  mineral  processing,  and  the  overall  uranium  extraction 
of  90.73  percent  could  be  attained. 
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SOURCE:  Beijing  YOU  KUANG  YE  [URANIUM  MINING  AND  METALLURGY]  in  Chinese  Vol  4, 
No  3,  Aug  85  PP  19-25 

TEXT  OF  ENGLISH  ABSTRACT:  Based  on  the  research  work  of  dissolution  rate  of 
U0„  in  various  systems  of  Mn02-U02-H2S0^-H20,  Mn02-Fe(III)-U02“H2S0^-H20  and 

U0__Fe ( III ) -H2SO11-H2O  with  Mn02  as  the  oxidant,  the  relationship  is  explored 
between  the  oxidation-dissolution  rate  of  UOp  and  the  dissolution  rate  of 
ferric  active  complex  formed  on  the  MnOp  surface,  and  the  electron  pair 
quantity  of  Fe(III)/Fe(II)  in  the  solution  as  well. 

It  is  believed  that  dissolution  rate  of  UOg  is  dominated  by  that  of  ferric 
active  complex  on  the  MnO_  surface.  In  the  fixed  conditions  of  leaching 
temperature  and  acidity  with  high  total  ferroferric  concentration,  the 
dissolution  rate  of  primary  uranium  mineral  is  mainly  dependent  on  the  Fe(III) 
concentration  in  leaching  solution  and  bears  no  important  relation  redox 
potential.  If  the  total  ferroferric  concentration  in  leaching  solution  is 
over  5g/1  (Fe(III)/Fe(II)=1 ) ,  the  redox  potential  as  low  as  -400  mV  or  even 
somewhat  lower  is  still  adequate  to  be  the  oxidative  requirement  for  leaching 
of  uranium  ore.  However,  if  the  total  ferroferric  concentration  in  leaching 
solution  is  relatively  low,  the  dissolution  rate  of  primary  uranium  mineral  is 
mainly  dependent  upon  the  ratio  of  Fe(III)/Fe(II) ,  i.e.,  the  redox  potential 

of  the  leaching  solution.  For  example,  when  the  total  ferroferric 
concentration  is  1g/1,  the  primary  uranium  mineral  can  be  well  dissolved  only 
if  the  Fe  (III)/Fe(II)  in  the  leaching  solution  is  maintained  at  >  10,  i.e., 

the  redox  potential  is  all  along  maintained  at  <  -400  mV. 
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TEXT  OF  ENGLISH  ABSTRACT:  In  this  paper  the  effects  of  the  main  influential 
technological  parameters  on  uranium  adsorption  process  in  the  multi-stage 
agitating  tanks  are  described.  These  parameters  are  resin-saturation,  resin- 
residence  time,  residence  time  of  processed  solution  and  uranium  concentration 
of  the  feed  solution.  The  relationship  among  these  parameters  is  presented  in 
the  expression  of  material  balance  equation. 

1.  Effects  of  uranium  concentrations  of  the  feed  solutions.  The  test  data 
indicated  that  the  variant  uranium  concentrations  of  the  feed  solutions,  i.e., 
at  the  range  of  0. 12~0.50g/1U,  did  not  make  a  numerical  difference  to  the 
number  of  adsorption  stages  required,  however,  the  specific  equipment  capacity *Y 
and  the  specific  adsorption  rate  6)  increased  with  the  increase  of  the  uranium 
concentration  of  feed  solutions  in  question. 

2.  Effects  of  residence  time  of  processed  solution.  At  the  range  of  residence 
time  from  5  to  30  minutes,  the  increase  of  residence  time  had  almost  no  effect 
on  ca  and  the  number  of  adsorption  stages,  but  the  T  increased  with  decrease 
of  the  residence  time  of  the  processed  solution. 

3.  Effects  of  resin-residence  time.  As  the  resin-residence  time  was  set  at 
40,  70,  100  and  200  minutes,  the  longer  the  resin-residence  time,  the  less  the 
adsorption  stages  needed  and  the  (£  ,  however,  the  higher  the  T . 

4.  Effects  of  resin— saturation .  The  tests  showed  that  the  higher  the  resin 
saturation,  the  more  the  adsorption  stages  required  and  the  lower  the  T  and 
the  <»  . 
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TEXT  OF  ENGLISH  ABSTRACT:  There  has  been  no  design  principle  so  far  available 
or  found  in  the  literature  for  engineering  design.  Thus,  a  practical  and 
useful  design  principle  for  the  mixed-phase  ports  of  box-type  mixer-settlers 
is  urgently  needed  for  extraction  processes.  In  this  paper,  hydraulic 
behaviors  taking  place  in  some  mixed-phase  ports  on  condition  that 
hydrostatic  balance  had  been  established  between  the  mixer  and  the  settler 
were  analyzed,  and  also  the  rule  of  distribution  of  hydrostatic  pressure 
difference  along  the  vertical  direction  in  the  mixed-phase  ports  was  revealed. 
The  causes  of  recycle  in  mixed-phase  ports  were  explained,  and  the  equations 
for  design  of  calculating  size  were  derived  for  various  cases.  It  was 
quantitatively  shown  that  the  mixed-phase  ports  should  be  designed  for  a  long 
horizontal  opening,  when  pressure  drop  might  be  negligible  and  there  was  not 
any  recycle. 
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